The Project Gutenberg EBook of Leçons sur l'intégration des Équations 
Différentielles aux Dérivées Partielles, by Vito Volterra 

This eBook is for the use of anyone anywhere at no cost and with 
almost no restrictions whatsoever. You may copy it , give it away or 
re-use it under the terms of the Project Gutenberg License included 
with this eBook or online at www.gutenberg.org 



Title: Leçons sur l'intégration des Équations Différentielles aux Dérivées Partielles 
Professées a Stockholm sur l'Invitation de S. M. le Roi de Suéde 

Author: Vito Volterra 

Release Date: August 24, 2009 [EBook #29783] 

Language : French 

Character set encoding: ISO-8859-1 

*** START OF THIS PROJECT GUTENBERG EBOOK ÉQUATIONS DIFFÉRENTIELLES *** 



Produced by Laura Wisewell, Andrew D. Hwang, and the Online 
Distributed Proofreading Team at http://www.pgdp.net (The 
original copy of this book was generously made available 
for scanning by the Department of Mathematics at the 
University of Glasgow.) 



Notes sur la transcription 

Ce livre a été préparé à l'aide d'images fournies par la Département 
de Mathématiques, Université de Glasgow. 

Des modifications mineures ont été apportées à la présentation, 
l'orthographe, la ponctuation et aux notations mathématiques. Les 
illustrations ont été recrées. 

Ce fichier est optimisé pour imprimer, mais peut être aisément 
reformaté pour être lu sur un écran. Veuillez consulter le préambule 
du fichier DT^X source pour les instructions. 



LEÇONS 



SUR 



L'INTEGRATION 



DES 



EQUATIONS DIFFERENTIELLES 

AUX 

DÉRIVÉES PARTIELLES 

PROFESSÉES À STOCKHOLM 

(FÉVRIER, MARS 1906) 

SUR L'INVITATION DE S. M. LE ROI DE SUÈDE 



PAR 

M. V. VOLTERRA 

SÉNATEUR DU ROYAUME D'ITALIE 
PROFESSEUR DE PHYSIQUE MATHÉMATIQUE À L'UNIVERSITÉ DE ROME 



UPSAL 1906 

IMPRIMERIE ALMQVIST & WIKSELL 



îere 



leçon. 



Introduction. 

1. Le cours que je ferai se rapportera à quelques points de la théorie des 
équations différentielles de la physique mathématique. On sait que la physique 
mathématique traverse une période de crise. On abandonne certaines idées pour 
en suivre de nouvelles. Tous ceux, par exemple, qui ont lu les éloquentes pages 
que M. Poincaré a consacré à cette question et ceux, qui ont pris connaissance 
de l'état actuel de la science dans le bel ouvrage de M. Picard, sont renseignés 
d'une manière fort claire là-dessus. Mais, même si certains concepts que nous 
avons maintenant sur la nature des phénomènes naturels et quelques principes 
fondamentaux devaient être ébranlés par de nouveaux faits et de nouvelles décou- 
vertes, une partie de la physique mathématique a bien des chances de se sauver 
du naufrage. Elle représente en effet, peut-être d'une manière grossière, mais cer- 
tainement d'une manière très-simple, une grande partie des faits naturels connus, 
les relie ensemble et a une utilité pratique hors de toute discussion. 

L'histoire des sciences nous offre l'exemple de théories analytiques de certains 
phénomènes qui ont été créées sous l'influence de certains principes et qui ont 
résisté à la chute de ces principes. Pour n'en citer qu'un seul, parmi la foule de 
ceux qui se présentent, il suffit de rappeler la théorie des instruments optiques qui 
se conserve dans ses lignes générales, tandis que les principes de l'optique ont subi 
tant d'évolutions. 

2. Les théories de la propagation de la chaleur, de l'hydrodynamique, de l'élas- 
ticité, des forces newtoniennes et de l'électromagnétisme peuvent être aujourd'hui 
traitées sous un point de vue commun, de sorte qu'elles peuvent constituer un seul 
chapitre d'analyse. On peut systématiser les méthodes qu'on emploie et classifier 
tous les faits qui s'y rapportent par la classification des équations différentielles 
dont ils dépendent. Cela conduit à envisager trois types d'équations différentielles 
qu'on a l'habitude d'appeler elliptique, hyperbolique et parabolique et des types 
mixtes. 

Lorsqu'on étudie les choses sous cet aspect une seule formule est capable 
de nous révéler sous une forme unique des propriétés qui se rapportent à des 
phénomènes en apparence divers entre eux. Quelquefois l'analogie est une simple 
analogie analytique, quelquefois elle va bien au delà. 

Supposons maintenant pour un instant que la base des faits physiques sur 
laquelle pose l'édifice analytique vienne manquer, ou qu'on la néglige. Cet édifice 
est tellement solide et utile par lui-même qu'il continuerait à subsister comme un 
des plus beaux chapitres de l'analyse. 

3. Le point de départ de toutes les considérations dont nous venons de parler 
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est d'envisager un ensemble continu ou un domaine à une, deux ou trois dimen- 
sions. A chaque point du domaine correspond une quantité scalaire ou un vecteur 
ou plusieurs quantités scalaires et plusieurs vecteurs liés par les équations dif- 
férentielles. Ces quantités sont quelquefois constantes par rapport au temps et 
quelquefois variables. Dans ce cas on a en général un grand avantage en consid- 
érant le temps comme une nouvelle coordonnée. 

Tout le monde sait que les idées des physiciens ont toujours oscillé entre le 
concept d'un milieu continu, siège de tous les phénomènes par lequel on tâche 
de supprimer toute action à distance, et le concept fondé sur l'hypothèse des 
molécules séparées et des actions à distance. 

Il ne faut pas croire que nos considérations soient liées forcément au pre- 
mier concept. Elles correspondent aussi à l'autre. Il suffit pour cela de rappeler 
que Cauchy, Poisson, Fourier, Laplace, qui suivaient les idées, qu'on appelle main- 
tenant de la mécanique physique, c'est à dire au fond le second système, ont été les 
premiers à découvrir les équations différentielles qui forment la base des théories 
analytiques. Il leur a fallu, pour y parvenir, faire un passage à la limite qui les 
a amené du discontinu au continu. Mais une fois cette limite franchie, les deux 
conceptions au point de vue analytique se mêlent dans la plupart des cas. 

4. Je n'aurai pas le temps de traiter d'une manière complète le chapitre 
d'analyse dont je viens de parler. Je toucherai seulement à quelques points, qui 
me semblent avoir un certain intérêt et que je crois nouveaux. Le rôle que certaines 
solutions polydromes jouent dans les différents cas, voilà un point que je tâcherai 
d'examiner avec quelque détail. Nous envisagerons ces solutions dans les différents 
types d'équations et pour nous familiariser avec elles nous étudierons d'abord les 
solutions polydromes dans le cas de l'équilibre élastique. Elles nous conduisent à 
des résultats pratiques et frappants qu'on peut montrer par des modèles matériels 
qui nous dévoilent leur vrai caractère et leur importance. Le cas d'équilibre des 
corps élastiques à connexion multiple, non assujettis à des forces externes, est 
en dépendance étroite avec les solutions polydromes, et nous offre de nouveaux 
problèmes de la théorie de l'élasticité. 

5. La théorie des fonctions est liée aux problèmes de physique mathématique 
qui dépendent de l'équation de Laplace à deux variables. Si l'on envisage, par 
exemple, un voile liquide incompressible infiniment mince qui est en mouvement, 
tout théorème de la théorie des fonctions peut être interprété comme un théorème 
relatif au mouvement. Réciproquement toute propriété du mouvement peut être 
interprétée comme un théorème de la théorie des fonctions. On obtient cette re- 
lation par les fonctions conjuguées qu'il faut regarder comme la partie réelle et 
le coefficient de l'imaginaire dans la théorie des fonctions et en hydrodynamique 
comme le potentiel de vitesse et la fonction qui définit les lignes des courants. 

Mais cela est borné au cas de deux dimensions. Comment généraliser la chose 
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au cas d'un nombre quelconque de dimensions ? Je montrerai que la théorie des 
fonctions conjuguées peut s'étendre au cas général par l'introduction d'un nouveau 
concept analytique sur lequel je reviendrai tout à l'heure. Cela nous amènera à 
exposer certaines vues nouvelles sur la théorie générale des fonctions. 

6. Par rapport aux équations du type hyperbolique j'exposerai d'abord, sans 
entrer dans les détails, quelques résultats que j'ai trouvés et publiés il y a quelque 
temps et qui plus récemment ont été étendus et complétés par d'autres géomètres. 

Ensuite je tâcherai de montrer le rôle bien singulier que joue le principe des 
images. La mémorable méthode de Lord Kelvin peut être employée quelquefois 
même dans ce cas et amène à des résultats plus simples que dans celui des équa- 
tions du type elliptique. L'influence de la réalité des caractéristiques se révèle par 
là d'une manière frappante. 

Il y a un mémoire très-profond de Weierstrass sur l'intégration des équations 
différentielles linéaires aux dérivées partielles où il expose une méthode par laquelle 
il intègre quelques équations du type hyperbolique. M me Kowalewski a appliqué 
dans un travail sur l'optique la méthode de Weierstrass, mais, sans s'en douter, elle 
avait à faire avec des fonctions polydromes, c'est pourquoi la méthode ne l'a pas 
amenée au résultat qu'elle cherchait. J'entrerai dans la discussion de cette question 
et du rôle de la solution polydrome. Je n'ai connu aucun géomètre, après M me 
Kowalewski, qui ait employé la méthode de Weierstrass et elle paraît comme une 
méthode à part qui n'est pas reliée aux autres. Je serai heureux de montrer qu'elle 
se rattache d'une manière très-simple à la méthode de Kirchhoff. Or, puisque les 
méthodes de Kirchhoff de Green et de Riemann ont des rapports étroits entre 
elles, toutes ces différentes méthodes restent reliées ensemble. 

7. Les dernières leçons seront consacrées aux équations différentielles du type 
parabolique. Leur étude n'est pas si avancée que celle des équations des autres 
types. Au premier abord il ne semble pas que la méthode des caractéristiques, qui 
a donné les résultats les plus généraux dans le cas hyperbolique, puisse conduire 
à embrasser tous les résultats connus lorsqu'on l'applique aux équations du type 
parabolique. Nous montrerons où est la difficulté. On a toujours conçu la méthode 
de Riemann comme bornée au domaine des variables réelles. Au contraire pour 
approfondir avec succès la question des équations du type parabolique, il faut 
commencer par étendre la méthode aux variables complexes et ensuite l'employer 
ainsi généralisée. 

Comme application des formules qu'on a dans les cas des équations du type 
parabolique j'exposerai la solution d'un problème qui se rapporte aux oscillations 
des fluides incompressibles. J'aurai ainsi l'occasion de dire quelques mots par 
rapport à cette question générale qui d'un côté est liée aux équations du type 
elliptique et d'un autre côté a bien des rapports avec les problèmes des ondes 
qu'on étudie par des équations du type hyperbolique. Je donnerai à ce propos une 
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formule analogue à celle ordinaire de Green qui contient une fonction analogue à 
celle qu'on appelle fonction de Green. 

8. Voilà le programme du cours. Cependant je n'ai pas encore parlé d'une 
manière explicite d'une question dont j'aurai l'occasion de traiter avec quelque dé- 
tail. Nous avons vu que dans nos considérations nous prenons comme fondement 
un continu à un certain nombre de dimensions, dans lequel on envisage des quan- 
tités scalaires ou des vecteurs. Supposons que nous ayons une quantité scalaire 
qui soit une fonction des points d'un domaine à trois dimensions, et pour fixer 
les idées supposons que ce soit la température des points d'une surface. Aura-t-on 
examiné d'une manière complète la question en envisageant cette quantité scalaire 
comme une fonction de trois variables, c'est à dire des coordonnées des différents 
points de la surface et du temps ? Il est facile de se convaincre que, si nous pou- 
vons changer arbitrairement la température au contour du corps, la température 
dépendra aussi de toutes les valeurs de la fonction qui exprime la température 
au contour. De même le potentiel d'un corps déformable dépendra de la forme 
du corps. On est par là amené d'une manière fort naturelle, et l'on peut même 
dire qu'on est forcé d'envisager non seulement les fonctions qui dépendent d'un 
certain nombre de variables, mais aussi celles qui dépendent de la forme de cer- 
taines lignes et de certaines surfaces et de toutes les valeurs de certaines fonctions. 
Nous aurons l'occasion de parler de ce concept et de développer quelques idées 
sur ces fonctions et sur leur inversion, lorsque nous entrerons dans le sujet dont 
nous avons touché au S 5. 
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1. Envisageons dans un plan les axes coordonnés x, y et la fonction 9 

arc tg — . 
x 
9 est l'angle que le rayon vecteur fait avec 

l'axe x. Partons d'un point M situé sur l'axe x, 
et, après avoir parcouru un cycle M NM' autour 
de l'origine, revenons au point de départ M\. 

En prenant les valeurs de 9 qui se suivent avec 
continuité la valeur 9\ qu'on trouve en M x après 
avoir parcouru le cycle est égale à la valeur ini- 
tiale augmentée de 2tt. Voilà un exemple très-simple 
d'une fonction polydrome. 

Les dérivées partielles — = -, — = 

' ' y oy 




ainsi que les dérivées 



dx x 2 + y* oy x* + y* 

successives sont finies continues et monodromes excepté à l'origine. C'est pourquoi 

f)9 f)9 
dans un domaine S qui ne comprend pas l'origine les dérivées — — , — - seront 

ox oy 

régulières ainsi que leurs dérivées successives. 

2. Rappelons un théorème de calcul intégral 
qu'on appelle théorème de Gauss. 

Soit S un domaine à trois dimensions X , Y , 
Z , trois fonctions monodromes finies et continues 
qui ont des dérivées monodromes finies et continues. 
On a 



dX 
dx 



OY 
dy 



dZ 

dz 



dS 



I (X cos nx + Y cos ny + Z cos nz) do 

•'a 




o étant le contour du domaine S, n la normale à a dirigée vers l'extérieur du 
domaine S. 

Rappelons aussi le théorème de Stokes. 

a étant une surface ayant pour contour S et X , Y , Z étant des fonctions 
monodromes finies et continues dont les dérivées sont aussi monodromes finies et 
continues on a 

' (dZ dY\ fdX 3Z\ {ÔY dX\ 

-, — — cos nx + ^r — cos ny + ( — ^— ) cos nz > do 

a { \ oy oz J \oz ox J 



dx dy 

± (Xdx + Ydy + Zdz) 
Js 



n étant la normale à a. On peut faire l'intégration sur le contour s dans une 
direction telle que l'intégrale du second membre soit précédée du signe +. 

3. Un domaine à trois dimensions est acyclique ou à connexion simple si 
toute ligne fermée que l'on peut tirer à l'intérieur peut se réduire infiniment petite 
par une déformation continue sans sortir du domaine. 

Si cette condition n'est pas vérifiée, le domaine est cyclique ou à connexion 
multiple. Supposons que par une coupure transversale un espace cyclique devienne 
acyclique, on dit alors que cette connexion est double. Si l'espace cyclique devient 
acyclique lorsqu'on fait deux coupures, on dit que la connexion est triple, et ainsi 
de suite. Dans tout domaine acyclique une ligne fermée peut être regardée comme 
le contour d'une surface située à l'intérieur du domaine. 

4. Dans la théorie de l'élasticité on peut envisager un milieu continu et 
l'on peut supposer qu'il existe un état naturel où il n'y a aucune tension interne. 
Déplaçons infiniment peu les points du milieu, et soient u, v, w les composantes 
des déplacements, alors la déformation du milieu est définie par les quantités que 
les Anglais appellent strains 

dw 

dz 

du dv 

732 = 723 = TT + TT, 713 = 731 = Tj 1" ~^~ , 721 = 7l2 = TT" + ~^~ ■ 

dz dy dx dz dy dx 

Il peut y avoir déplacement sans déformation lorsque ces quantités sont nulles. 

Le potentiel des forces élastiques relatif à chaque élément dS du milieu dé- 
formé s'exprime par 

^(7ll, 722, 733, 723, 731, 712) dS 

où F est un polynôme homogène de 2 d degré par rapport à 7n, 722, 733, 723, 731, 
712 qui est toujours négatif et s'annule seulement lorsque toutes les quantités j rs 
sont nulles. 

La tension dans chaque point est caractérisée par les quantités 

dF dF dF 

tu — — , £22 — — , £33 







du 




dv 






7n = 


dx 


722 = 


dy 


733 : 


dv 


dw 






dw 


du 



O711 C>722 0733 

dF dF _ dF 

£32 — ^23 — ô , ^13 — £31 — 7^ , £21 — £l2 



C>732 C731 C>7i2 

que les Anglais appellent stress. 

La tension unitaire qui s'exerce sur chaque élément d'aire do du milieu a 



pour composantes 

T x = t u cos nx + ti2 cos ny + £ 13 cos nz 
T y = £21 cos nx + £22 cos ny + £23 cos nz 
T* = £31 cos nx + £32 cos ny + £33 cos nz 

n étant la normale de l'élément do. 

Si X, Y, Z sont les composantes des forces extérieures unitaires qui s'exercent 
sur chaque élément dS du corps élastique, et T x , T y , T z les composantes de la 
tension unitaire qui s'exerce sur chaque élément du contour, les équations de 
l'équilibre élastique sont 



(1; 



X 

Y 
Z 



dtn , dt 
dx 
dt 21 



12 



dt 



13 



dx 

dtu 
dx 



dy 
dt 22 

dy 
dt 32 

dy 



dz 

dt 23 

dz 

dt 33 

dz 



T x = £n cos nx + £12 cos ny 
(2) < T y = £ 2 i cos nx + £22 cos ny 
[ T z = £31 cos nx + £32 cos ny 



où n est la normale au contour dirigée vers l'intérieur du corps. 
Lorsque le corps est isotrope et homogène on a 



£13 cos nz 
£23 cos nz 
£33 cos nz 







F = \LB^ 


+ X^, 




en ayant posé 










= 7n 


+ 722 + 733, tp = 7ll + 


7 2 2 2 + 7 3 2 3 + 


m, 


c'est pourquoi 


£ll 
£22 
£33 


= L9 + 2X 7ll 
= L6 + 2X722 
= L6 + 2X733 


£23 = 

£31 = 

£12 = 


-- Kl2Z 

-- Kj 31 
'- K712 



il 



7? 2 ); 



L et X sont des quantités constantes et négatives, 9 est la dilatation cubique. 
Les équations (1) deviennent alors 



(!') 



X = XA 2 n + (L + X) 



y 



XA 2 t> 



(L + X) 



Z = KA 2 w + (L + X) 



90 

de 

dy 

de_ 

dz 



le symbole A 2 représentant l'opération 

d 2 d 2 d 2 



2 /-) *-z 



dx 2 dy 2 dz 

6. Ayant rappelé ces notions fondamentales nous voulons montrer que la 
théorie de l'élasticité pour les corps ayant une forme cyclique se présente d'une 
manière différente que pour les corps acycliques, et cela tient au rôle que les 
solutions polydromes jouent dans la question. 

On énonce en général le théorème qu'tm corps élastique qui n'est pas soumis 
à des actions externes est en équilibre, et on tire de là que les forces externes étant 
données la déformation du corps est connue. Or ces propositions ne sont exactes 
que si le corps a une connexion simple. 

Pour montrer qu'il y a des cas où elles ne se vérifient 

/-" -"--\ pas, envisageons un anneau. Retranchons une tranche ra- 

/ \ diale AA', BB' très mince, et soudons les deux bouts AA' 

[ \, \ et BB'. Le corps après la soudure sera soumis à des ten- 

[ =- ") ] sions internes et cependant il ne supportera aucune action 

\ V b/^-^Jb' extérieure. 

\ / On pourrait soupçonner que quelque discontinuité ou 

^\^_^^^ quelque singularité s'est produite dans la déformation à 

l'endroit où l'on a fait la soudure ; mais on peut démontrer 

qu'il n'y a aucune singularité de cette sorte, de manière que si l'on voulait trouver 

l'endroit où l'on a fait la soudure par quelque particularité de la déformation il 

serait impossible de la retrouver. 

On peut se demander à quoi tient la contradiction entre les propositions 
précédentes et le résultat que nous venons d'examiner. 

Le nœud de la question est dans l'application du théorème de Gauss. En 
effet pour démontrer les propositions que nous avons énoncées tout à l'heure on 
suppose que dans les équations (1) et (2) X = Y = Z = T x = T y = T z = 0, on 
multiplie les deux membres des équations (1) respectivement par u, v, w, qu'on 
somme, qu'on intègre et transforme par le théorème de Gauss. On tire de là que 

FdS = 0, 

s 

S étant l'espace occupé par le corps élastique ; d'où l'on déduit 

7ll = 722 = 733 = 723 = 731 = 7l2 = 

c'est à dire que si les forces extérieures sont nulles il n'y a pas de déformation. 
Or nous avons remarqué que pour appliquer le théorème de Gauss il faut que 



toutes les fonctions qu'on emploie soient monodromes. Montrons que dans le cas 
de l'anneau les composantes des déplacements u, v, w sont polydromes. 




En effet examinons les déplacements après la soudure. Les points de AA' 
ne se seront pas déplacés tandis que les points de BB' se seront déplacés de la 
largeur de la fissure qu'on a faite. C'est pourquoi si nous partons d'un point M 
de AA', parcourons un cycle à l'intérieur de l'anneau et prenons les valeurs des 
déplacements qui se suivent avec continuité, nous revenons au point de départ de 
l'autre côté M x de la coupure avec des valeurs des déplacements différents des 
valeurs qu'on avait au point de départ. C'est un cas de polydromie analogue à 
celui que nous avons envisagé dans le premier §. 

On tire de là deux conséquences : 

La première est que pour rendre exactes les deux propositions qu'on a énon- 
cées il faut ajouter la condition que les déplacements soient monodromes. 

La seconde conséquence est qu'il faut reposer la question : 

Quand est-ce que les déplacements peuvent être polydromes ? 

7. Dans les considérations que nous venons de faire nous nous sommes laissés 
conduire par l'intuition ; mais il serait dangereux de continuer de la sorte. En effet 
par intuition on serait amené à penser que le même résultat qu'on a obtenu dans le 
cas de l'anneau on pourrait l'obtenir dans le cas d'un corps acyclique quelconque 
en y retranchant une tranche très-mince et en soudant les faces de la coupure. Or 
dans ce cas, par cette opération, la déformation du corps ne serait pas régulière, 
c'est pourquoi tout se présenterait d'une manière différente, que dans le cas de 
l'anneau. 

Nous allons démontrer que si l'on suppose que la déformation est régulière, 
c'est à dire que si 711,722,733,723,731,712 sont des fonctions finies monodromes 
et continues et que leurs dérivées sont aussi finies continues et monodromes, la 
polydromie des déplacements ne peut se présenter que si le corps a une forme 
cyclique. 

Il suffit de trouver des formules par lesquelles on puisse calculer u, v, w en 
connaissant les quantités j rs . 

Soit 5 une ligne, et désignons par A et Ai ses points extrêmes. Soient 7rs 
et 7rs les valeurs de 7 rs aux points A et Ai, x y zq les coordonnées de A , xi yi zi 



(I) 



(I") 



les coordonnées de A-y. Les valeurs de u, v, w au point Ai seront données par les 
formules 

1^(121 + r )(yi - yo) + 2 1 

dx 



Ui = u 
+ 



lOl = W 

+ 



;(7l? 



Tfii + (Vi-v)-£— + (zi-z) 



(yi - y) 
yi-y 



dy 

<9y <9x 
^72i , <973i 



?o)(*i -*o) 

dx 
ds 



dz 

Z\ — z 



«9721 , 9731 07: 



23 



(9.3 



«97 



23 



cb 



<9y 



<9x - 



+ (zi - z) 



<9y 

^713 



<9x 

^733 
<9x 



(!') f 1 = w + 



,(°) 



(732 ; + Po)(*i - zo) + -(7S - r )(xi - x ) 



Z\ — z 



«9732 , «97i2 (97 



31 



dx 



dz 



722 + (21 - z) — h (xi - x) 



c'y 

•9722 



+ (xi - x) 



11 



(972i _ (97: 
<9x <9y 



( Zl - z) 



dz 

•9723 «9733 



<9z 



<9y 



dx 



dy 
ds 



«9732 , (9712 ^7; 



31 



dx 



dz 



dy 



1 



(7Ï? + 9o)(»i - &o) + ^(7 2 3 J - Po)(î/i - yo) 



(o) 



2 

(xi — x) 



«9731 $7 
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yi-y\ (dl\z , «9723 <97 



<9x dz J \ 2 

^713 , «9723 ^712 
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c'y 



<9x 



dz 

^733 



733 + (*i-*)^ + (yi-t/) Ô2/ 



dy 

+ (2/1 - 2/) 

cl 2; 



<9x 

(9732 

«9y 



(92 
(9722 



ds 



ds 



dy 
ds 

dz 
ds 



dx 
ds 



dz 
ds 



d,x 
ds 
dy_ 

ds 



ds 



ds 



où u , v , w sont les valeurs de u, v, w au point A et p , q , r sont égales aux 
valeurs des différences 

dv dw dw du du dv 
dz dy ' dx dz' dy dx 

au point A . 

On le vérifie aisément en faisant les dérivées de ui, i>i, W\ par rapport à Xi, 

Vi, zi- 

Il faut maintenant chercher les conditions pour la monodromie. 

Supposons que la ligne s se réduise à un cycle fermé. Si le corps élastique oc- 
cupe une région acyclique, on pourra regarder 5 comme le contour d'une surface a 



située à l'intérieur du corps ; c'est pourquoi les intégrales qui paraissent dans les 
formules précédentes pourront se transformer par le théorème de Stokes et elles 
s'écriront 



2/1-2/ 



B 



Z\ — z 



C 



( Zl -z)E+^-^B 



{yi -y)E+^-^C 



t w , yi ~ y a 

(zi - z)F + -^—A 

X\ — X 



Z\ — z 



C 



[x\ — x)t H — G 



où n est la normale à a et 



A 

B 

C 



d /<973i «9712 «97: 



23 



dx \ dy 
d /Ô712 



dy \ dz 



dz dx 

<9723 _ ^731 

dx dy 
dj3i <97i2 



d 2 ju 
dzdy 

d 2 j22 
dxdz 

9 2 733 



cos ny - 
cos ny - 
cos ny 

E = 
F = 
G = 



(yi-y)G+^-^A 



{ Xl -x)G+^—^B 



2 
2/1-2/ 



B 



cos nz > du 
cos nz > da 
cos nz > da 



«9 2 732 «9 2 722 <9 2 733 



dydz dz 2 



dy 2 



<9 2 73i (9 2 7 33 (9 2 7n 



dzdx dx 2 dz 2 
<9 2 7i2 (9 2 7n <9 2 722 



dz \ dx ' dy dz ) ~ dydx dxdy dy 2 dx 2 

Mais par un théorème de De Saint Venant on a 

A= B =C = E = F =G=0 
c'est pourquoi les intégrales des formules (I)(I')(I") s'annuleront, et on aura 

U\ = Uq V\ = Vq W\ = Wq 

ce qui prouve que si la forme du corps est acyclique et la déformation régulière, 
les déplacements sont monodromes. Mais si le corps est cyclique, le cycle fermé s 
n'est pas en général le contour d'une surface a appartenant à la région occupée 
par le corps, et par suite on ne peut pas en général appliquer le théorème de Stokes 
et enfin les déplacements peuvent être polydromes. 

8. On peut résumer les résultats qu'on vient d'obtenir en énonçant les propo- 
sitions suivantes. Un corps élastique ayant une forme acyclique et une déformation 
régulière peut être amené à l'état naturel par des déplacements finis continus et 
monodromes. 

Si le corps élastique a une forme cyclique, cela n'est pas toujours possible, et 
pour l'amener à l'état naturel il sera nécessaire quelquefois de faire des coupures 
et de retrancher des parties du corps. 

Un corps élastique fini ayant une forme acyclique et n'étant pas sujet à 
des forces extérieures est à l'état naturel si la déformation ne peut cesser d'être 



régulière, mais si le corps est cyclique, le corps peut être à l'état de tension, même 
si n'existent pas des actions extérieures. 

Les forces extérieures étant connues, la déformation régulière d'un corps acy- 
clique est connue, mais celle d'un corps cyclique n'est pas connue. Elle est connue 
seulement dans le cas où l'on sait d'avance que le corps élastique peut être amené 
à l'état naturel par des déplacements monodromes. 



neme 



leçon. 



up - 


-u a 


= l + ry — qz 


V/3- 


-V a 


= m + pz — rx 


■-V/3 ~ 


-w a 


= n + qx — py 



9. Lorsque le corps est cyclique, on peut le réduire acyclique en faisant des 
coupures. Les formules (I)(I')(I") donnent aisément la loi des discontinuités des 
déplacements le long des coupures. Les valeurs des déplacements d'un côté d'une 
coupure étant u a , v a , w a et celles de l'autre côté étant up, vp, wp on a 



(3) 



l, m, n, p, q, r, étant des quantités constantes pour chaque coupure. 

A chaque coupure correspondent donc 6 valeurs constantes qu'on peut appeler 
les constantes de la coupure, et l'on a le théorème : 

Un corps élastique cyclique étant déformé régulièrement, la déformation est 
déterminée par les forces extérieures et les constantes de chaque coupure. 

On voit par là que les problèmes de l'élasticité pour les corps cycliques se 
présentent d'une manière fort différente que pour les corps acycliques. 

10. Les formules (3) montrent très-facilement la signification mécanique des 
constantes de chaque coupure. En effet supposons que nous fassions effectivement 
les coupures dans le corps cyclique déformé de manière qu'il puisse reprendre l'état 
naturel; et s'il est nécessaire retranchons les parties surabondantes du corps. 

Alors les formules (3) nous montrent que les molécules situées des deux côtés 
de chaque coupure et qui adhéraient auparavant, ont été séparées par un déplace- 
ment relatif qui résulte d'une translation et d'une rotation. Cette translation et 
cette rotation sont égales pour tous les couples de molécules qui adhéraient le long 
de la coupure. 

Si l'on prend pour centre de réduction l'origine, les constantes de la coupure 
sont les composantes de la translation et de la rotation par rapport aux axes coor- 
donnés. 

Réciproquement pour déformer le corps cyclique on pourra faire les coupures 
qui le rendent acyclique. Ensuite, l'on pourra déplacer les deux faces de chaque 
coupure de manière que le déplacement relatif de tous les couples de molécules qui 



se trouvent vis à vis soit résultant de la même rotation et de la même translation. 
Enfin on pourra souder entre elles les faces en retranchant ou en ajoutant la 
matière qui sera nécessaire. 

Nous appelons cette opération une distorsion du corps et les 6 constantes de 
chaque coupure les caractéristiques de chaque distorsion. 

Les formules (I)(I')(I") conduisent aussi immédiatement à un autre résultat. 

Prenons deux coupures qui peuvent se réduire l'une à l'autre par une trans- 
formation continue, alors les caractéristiques des deux coupures sont égales entre 
elles. 

On voit par là que les caractéristiques d'une distorsion ne dépendent pas 
des coupures que l'on a faites, mais de la déformation du corps et de sa nature 
géométrique. Le nombre des distorsions indépendantes d'un corps est égal à l'ordre 
de la connexion diminué d'une unité. 

11. Une question qui se pose naturellement est la suivante : Peut-on choisir 
toutes les caractéristiques d'une manière arbitraire ? 

On démontre en général que les caractéristiques des distorsions étant données 
d'une manière arbitraire, le problème de la déformation du corps peut être ramené 
à un problème ordinaire de l'équilibre où le corps n'est assujetti à aucune distorsion 
mais seulement à des forces données. 

Nous laisserons de côté ce théorème, général, et nous envisagerons un cas 
particulier. Soit a une aire finie du plan xz qui ne rencontre pas l'axe z. Faisons 
tourner le plan autour de cet axe et en même temps déformons et déplaçons l'aire a 
dans le plan de manière qu'elle ne rencontre jamais l'axe z et qu'elle revienne à la 
forme et à la position initiale après avoir fait un tour. 

L'aire a décrit un volume à connexion double qu'on peut supposer rempli 
d'une substance élastique isotrope. Z, m, n, p, q, r, étant des constantes, posons 

1 \n ... , ...a ....... V , 1 ( ... .... rK _\,__, 2 , 2n 



u= — Ul-qz + ry) arc tg - + - ( -m - pz - L + 2K x j lo &( x +V ) 

lf y 1 / rK \ 

( 4 ) { v = —< (m- rx + pz) arc tg- + - [l- qz- y log(a; 2 + y 2 ) 

Ziï y x Z \ L + ZK J 

w = ^~ \ ( n ~ py + i x ) arc te ~ + o(p x + iv) l og(x 2 + y 2 ) 

ztï y x z 

On vérifie aisément que ces fonctions satisfont les équations (1') et en même temps 
à cause de la polydromie de la fonction arc tg - (Voir § 1) qu'elles sont polydromes 
et correspondent à la distorsion la plus générale du corps, l, m, n, p, q, r sont les 
caractéristiques de la distorsion. Il est facile de calculer les tensions qui agissent 
sur la surface libre du corps. 
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1. Le problème fondamental de la théorie des distorsions des corps élastiques 
ayant une forme cyclique est de déterminer la déformation du corps lorsque les 
caractéristiques des distorsions sont données. 

Nous allons calculer l'énergie d'un corps élastique qui n'est assujetti à aucune 
action externe et dont les distorsions sont connues. 
L'énergie est donnée par 

E = - I FdS = -- / (tn7n + £22722 + £33733 + £23723 + £31731 + £12712) dS 
Js z J s 

où S est le volume du corps. 

On peut maintenant transformer cette intégrale en employant le théorème de 
Gauss, mais il faut prendre garde que u, v, w sont des fonctions polydromes. C'est 
pourquoi l'application directe du théorème ne peut pas se faire comme nous avons 
déjà vu dans la leçon précédente. 

2. Pour pouvoir appliquer le théorème de Gauss il faut commencer par sec- 
tionner le volume S de sorte qu'il devienne acyclique, et il faut ensuite regarder 
les deux faces de chaque coupure comme faisant partie du contour du corps. 

On trouve alors, en désignant par <7i <t 2 ■ • ■ cr n les coupures, que la transfor- 
mation de l'intégrale conduit à la formule suivante 

1 n r 

E= 2^J {X ^ ~ ^ + Yi{v « ~ ^ + Zt(w « ~ W ^ ))} dai 
= r^jU Xi da { + rrii Y t dai + n; / Z { da { 

I L J ai J <Ji J a. L 

+ pi (Y t z - Ziy) dai + q t (ZiX - X, t z) dai + n J {Xiy - Yix) dai \ 

où Ua -, Va, Wa ; uv, Vg\ Wg', sont les valeurs des composantes des déplacements 
des deux côtés de la coupure <7j : li,rrii,rii,'Pi,qi,r i sont les caractéristiques de la 
distorsion relatives à la coupure <7j ; et X t Yi Z { sont les composantes de la tension 
unitaire qui s'exerce sur chaque élément d'aire de la même coupure. 
Posons 

Li= X t dai Mi= Yi dai N t = Z t da { 

Pi= j (YiZ - Ziy) dai Qi = {Z { x - X { z) da { Ri = (X t y - Yix) da { 

-J ai ^ ai J ai 

1 n 

E = - y^X L ih + M i m i + N i n i + P iPi + Qili + R i r i)- 



on aura 



2 

i=l 
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Si l'on compose les tensions qui sont appliquées sur la coupure c^, en prenant 
pour centre de réduction l'origine, on trouve que Lj, Mi, Ni sont les composantes 
de la force résultante et Pi,Qi, Ri les composantes de la couple résultante. 

On peut écrire la formule précédente d'une manière plus simple en appelant 
Si s 2 • • • s en les quantités U, rrii, rii et E 1 , E 2 . . . E Qn les quantités correspondantes 
Li,Mi,Ni; elle devient en effet 

1 

Supposons que les quantités si s 2 ■ ■ ■ sq ti soient nulles excepté Sh = 1. Ap- 
pelons E ih les valeurs correspondantes des coefficients E im Lorsque s\ s 2 . . . s 6n 
sont quelconques on aura 

(m 

Ei = 2_^ E ih Sh 
1 
c'est pourquoi 

_. 6n 6n 

e = ^y^y^ E ih sis h . 

i=l h=\ 

3. Nous allons démontrer maintenant une propriété fondamentale des co- 
efficients E^, et nous nous permettons à ce propos de faire quelques remarques 
générales. 

Green a démontré par l'application du théorème de Gauss un théorème fonda- 
mental dans la théorie du potentiel. Mais le théorème de Green n'est pas borné au 
cas du potentiel, il peut s'étendre à un nombre considérable de cas. J'ai démontré 
qu'il peut s'étendre à tous les problèmes qui dépendent du calcul des variations. 

Même le problème de l'élasticité, comme tout problème de mécanique, se 
rattache à une question de calcul de variation, c'est pourquoi on a le théorème 
analogue à celui de Green pour l'élasticité. Il a été donné pour la première fois 
par Betti qui en a fait des applications remarquables pour l'intégration des prob- 
lèmes de l'équilibre élastique. Cerruti, Boussinesq, Somigliania ont continué par 
le chemin qui avait été tracé par Betti. 

Or si les déplacements sont polydromes, puisque le théorème de Gauss ne 
peut plus s'appliquer, le théorème de Betti n'est plus applicable. Nous allons voir 
cependant qu'on peut reprendre, même dans ce cas, l'idée fondamentale de Green, 
et l'on est amené par là à une loi de réciprocité fort intéressante. 

Supposons que nous appliquons successivement au corps élastique deux dis- 
torsions caractérisées par les valeurs s[, s' 2 , . . . s' 6n et s", s' 2 ' . . . Sg„ des caractéris- 
tiques. Soient j' rs et j'^ les valeurs des j rs correspondant aux deux déformations, 
F' et F" les valeurs du potentiel élastique. 
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On aura, par un théorème bien connu d'algèbre, 



et si nous transformons ces intégrales, après avoir réduit acyclique le volume S par 
les coupures, en regardant toujours les faces des coupures comme faisant partie 
du contour du volume S on trouve 



6n 6n 



E £ « 



771/ // 



1 1 

E" et El étant les valeurs des quantités E t correspondant aux deux distorsions. 
On tire de là 

6n 6n 6n 6n 

i=l h=\ i=l h=l 

et puisque les quantités s^ et s" sont arbitraires, on doit avoir 

Ea = E hi . 

Cette propriété réciproque des coefficients E ih peut s'interpréter de plusieurs 
manières et conduit à des théorèmes mécaniques qu'on énonce facilement. 

On peut appeler l'ensemble de la force ayant pour composantes L*, Mj,JVj 
et de la couple ayant pour composantes Pi,Qi,Ri l'effort total appliqué à la 
coupure a« et l'on peut appeler toutes les quantités Ei, E 2 . . . E 6n les caractéris- 
tiques des efforts. 

On a alors que si deux systèmes de distorsions engendrent deux systèmes 
d'efforts, la somme des produits des caractéristiques des efforts du premier système 
multipliées par les caractéristiques des distorsions du second système est égale à la 
somme des produits des caractéristiques des efforts du second système multipliées 
par les caractéristiques des distorsions du premier système. 

Lorsque les quantités Si,s 2 ■ ■ ■ ,s en sont nulles excepté s^ = 1 la distorsion 
peut s'appeler une distorsion élémentaire d'ordre h. E^ est l'effort d'ordre i en- 
gendré par la distorsion élémentaire d'ordre h. Ehh est l'effort conjugué de la 
distorsion élémentaire. Le théorème précédent peut donc s'énoncer en disant que 
l'effort d'ordre i engendré par la distorsion élémentaire d'ordre h est égal à l'effort 
d'ordre h engendré par la distorsion élémentaire d'ordre i. 



iieme 



leçon. 



4. Nous avons démontré dans la leçon précédente que la même distorsion exé- 
cutée sur deux coupures qui peuvent se réduire l'une à l'autre par une déformation 
continue engendre la même déformation du corps. 
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Les deux coupures peuvent s'appeler à cause de cela des coupures équivalentes. 
Le théorème qu'on vient de rappeler peut être complété en démontrant que les 
efforts dans les coupures équivalentes sont égaux. Prenons en effet la partie du 
corps entre deux coupures équivalentes. Les tensions qui sollicitent les éléments 
des coupures doivent se faire équilibre. On tire de là tout de suite l'égalité des 
deux efforts. 

Donc les efforts dépendent, comme les distorsions, de la nature géométrique 
du corps et de la déformation. 

Le problème que l'on peut se proposer est d'étudier les propriétés des efforts 
et de les déterminer, les distorsions étant données. 

Nous allons voir que l'on peut obtenir des théorèmes généraux sans intégrer 
directement les équations de l'élasticité. 

5. Supposons que nous ayons un solide de révolution qui ait une connexion 
double ; il peut être engendré par la révolution d'une aire plane simplement con- 
nexe autour d'un axe du plan qui ne la rencontre pas, ou il peut être engendré 
par la rotation d'une aire doublement connexe limitée par l'axe. Supposons que le 
solide soit rempli d'une substance électrique dont la constitution soit symétrique 
par rapport à l'axe. 

L'énergie du système est donnée par 

6 6 

e = ^y^y^ E ih sis h 
i i 

où 

si = l, s 2 = m, s 3 = n, 5 4 =p, s 5 = q, s 6 = r. 

Cette expression peut se comparer à la force vive d'un liquide indéfini sans tour- 
billons à l'intérieur duquel se trouve un corps solide symétrique, et l'on peut 
répéter dans ce cas des raisonnements tout à fait semblables à ceux que l'on fait 
dans cette question d'hydrodynamique. On trouve par là, en vertu de la symétrie, 
qu'on peut toujours choisir l'origine en un point de l'axe de symétrie de manière 
que l'expression de E devienne 

E = l{E n (sî + si) + £33 s 2 3 + £44^ + 4) + E m si}. 

L'origine s'appelle le point central de l'axe de symétrie. 

Les quantités E ih sont nulles lorsque i et h ne sont pas égaux ; c'est pourquoi 
on a le théorème : 

Dans un corps symétrique qui a une connexion double chaque distorsion élé- 
mentaire engendre un seul effort qui est l'effort conjugué, le centre de réduction 
étant dans le point central de l'axe de symétrie. 
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Donc si la distorsion est due à une translation relative des molécules des 
deux faces de la coupure, l'effort total engendré est une force qui passe par le 
point central, et si la distorsion est une rotation l'effort total est une couple. 

Revenons à l'exemple primitif des distorsions. 
Prenons un anneau symétrique, retranchons une mince 
tranche radiale, et soudons les deux bouts. On 
serait facilement amené à croire que les faces de la 
soudure sont sollicitées par une traction, mais si nous 
réfléchissons que le déplacement relatif des molécules 
situées d'un côté de la coupure par rapport à celles 
situées de l'autre côté est dû à une rotation autour de 
l'axe de symétrie, on tire du théorème précédent que l'- 
effort est une couple ayant pour axe l'axe de symétrie. 

Mais si la résultante des tensions qui s'exercent sur les faces de la soudure est une 
couple, il faut qu'une partie de la soudure soit comprimée et qu'une partie soit 
tendue, et même que la somme des efforts de traction égalise la somme des efforts 
de compression. 

On arrive par là à un résultat qui est inattendu. 

Supposons maintenant que la même tranche re- 
tranchée ait une largeur uniforme. La distorsion con- 
siste alors dans une translation, et l'effort est une force 
normale à l'axe de symétrie et qui le rencontre au point 
central. Même dans ce cas une partie de la soudure est 
comprimée et une partie supporte une traction, mais si 
nous envisageons l'énergie du système, on a aisément 
que dans le cas de la fissure radiale la partie comprimée 
est la partie interne de la soudure et la partie externe 
supporte la traction, tandis qu'on a le contraire si la 
fissure est de largeur uniforme. 

Appliquons maintenant le théorème des coupures équivalentes. On en tire 
que si la fissure est radiale toute section de l'anneau supporte, d'une manière 
symétrique, la compression et la traction, mais dans le cas de la fissure uniforme 
dans la région A' de l'anneau opposée à la soudure A les choses sont renversées, et 
une section en A' supporte dans la partie interne la compression et dans la partie 
externe la traction. 

Voilà que par la considération de l'énergie du système déformé on tire un 
grand nombre de propriétés sans intégrer les équations de l'équilibre élastique. 

12. Il est intéressant cependant d'approfondir le cas de l'anneau symétrique 
isotrope en étudiant en détail les déformations produites par les distorsions. Il faut 
pour cela intégrer effectivement les équations de l'équilibre et il est nécessaire de 
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partir des formules (6) que nous avons données dans la précédente leçon. 

Nous envisageons un cylindre creux symétrique isotrope, nous prenons l'orig- 
ine dans le centre de symétrie, et faisons coïncider l'axe z avec l'axe de symétrie. 
Les déplacements donnés par les formules (4) de la leçon précédente correspondent 
à la distorsion la plus générale ayant pour caractéristiques l,m,n,p,q,r ; mais si 
nous calculons les tensions qui sollicitent le corps, nous voyons que pour obtenir 
ce cas d'équilibre il faut supposer que les surfaces latérales et les bases du cylindre 
soient sollicitées par des tensions. Or nous voulons envisager le cas où le cylindre 
est déformé par la distorsion et n'est pas assujetti à des actions externes. Il faut 
donc éliminer les tensions superficielles. 

Représentons par T l'ensemble de ces tensions. Supposons maintenant que le 
cylindre primitif n'ayant pas de distorsions soit soumis à l'action des tensions — T, 
c'est à dire des tensions égales et contraires aux tensions T. Soient u' v' w' les 
déplacements qu'on trouve. En prenant les différences 



u — u 



«• 



w 



on aura des déplacements qui correspondent à la distorsion du corps, mais en même 
temps les tensions au contour auront été éliminées. On voit par là que l'élimination 
des tensions au contour est un problème ordinaire de l'équilibre élastique qu'on 
peut tâcher de résoudre par les méthodes ordinaires. 

Dans les formules (4) nous avons 6 constantes arbitraires qui correspondent 
aux six distorsions élémentaires. On peut simplifier beaucoup le problème en en- 
visageant séparément chaque distorsion élémentaire. On voit par la symétrie que 
les distorsions d'ordre 1 et 2 se ramènent l'une à l'autre et de même celles d'ordre 4 
et 5. Il suffit donc d'envisager les distorsions d'ordre 

2, 3, 5, 6. 

Commençons par celles d'ordre 6 et 2. La distorsion 
d'ordre 6 correspond évidemment à la fissure radiale 
et celle d'ordre 2 à la fissure de largeur uniforme. 

Dans le cas de la distorsion d'ordre 6 on élim- 
ine très aisément, en vertu de la symétrie, les ten- 
sions latérales, et l'on trouve qu'après la distorsion 
le corps garde sa forme cylindrique régulière et les 
bases gardent leur distance initiale si on les com- 
prime dans la région (interne) que nous avons lais- 
sée en blanc et si on exerce une traction dans la 
région (externe) que nous avons hachée. Appelons T l'ensemble de ces tensions. Il 
n'y a plus de difficulté alors à déterminer la forme du corps lorsqu'il est complète- 
ment libre et qu'aucune force ne s'exerce sur les bases. En effet il suffit d'envisager 
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un cylindre à l'état naturel et de chercher la déformation lorsqu'on le sollicite par 
les tensions —T. La forme que le corps prendra sera celle que nous cherchons. 
Prenons une tranche radiale du cylindre : la partie en blanc supportera une trac- 
tion, la partie hachée une compression ; c'est pourquoi cette tranche subira une 
flexion. 






Chaque tranche du cylindre se déformant par une flexion pareille le cylindre 
prendra la forme indiquée par la figure voisine, c'est à dire que le bord interne 
se soulèvera, le bord externe s'abaissera tandis que la surface latérale externe 
deviendra concave, et la surface latérale interne convexe. 

J'ai voulu vérifier par l'expérience ce résultat, et j'ai pris un gros cylindre 
creux de caoutchouc, j'y ai fait une fissure radiale, et j'ai soudé les deux faces de 
la fissure. La forme prise par le corps a été celle que le calcul avait prévue, et l'on 
a pu même vérifier par des mesures que les résultats des calculs étaient exacts. 
On pouvait aussi remarquer que la partie interne de la soudure était comprimée 
et que la partie externe supportait un effort de traction. A cause de cet effort la 
soudure ne résistait pas, c'est pourquoi pour garder la forme j'en ai pris le modèle 
en plâtre que je vous montre. L'endroit de la coupure est parfaitement visible. 

Quelquefois il arrive que l'on veut rétrécir un tube, mais si l'on retranche une 
tranche radiale du tube et que l'on soude les deux faces, le tube prend sa forme 
cylindrique. Les forgerons et tous ceux qui ont essayé cette opération connaissent 
très-bien ce résultat que maintenant la théorie des solutions polydromes explique 
parfaitement. 

13. Quelle est maintenant la forme du cylindre lorsqu'on fait la fissure de 
largeur uniforme, c'est à dire lorsqu'on fait une distorsion d'ordre 2 ? 

On peut procéder tout à fait de la même manière que dans le cas précédent. 
On peut commencer par éliminer les tensions latérales et après les tensions sur 
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les bases. La forme prise par le corps est même plus compliquée que dans le cas 
précédent. Elle est représentée par la figure ci-jointe, et voici un plâtre qui est 
le modèle d'un cylindre de caoutchouc qui a supporté une distorsion d'ordre 2. 
Même dans ce cas l'expérience a vérifié complètement les résultats du calcul. 

Je n'entre pas dans les détails sur les 
distorsions d'ordre 5 et d'ordre 3. La pre- 
mière peut se déduire de celle d'ordre 2 
par des intégrations simples, et la dernière 
peut se calculer directement sans difficulté. 
Cette distorsion s'obtient en faisant glisser 
les deux faces de la coupure, l'une par rap- 
port à l'autre, dans la direction de l'axe du 
cylindre et en les soudant ensuite. 

Nous avons dit (§ 4) qu'on peut se 
proposer le problème de déterminer les ef- 
forts lorsqu'on connaît les distorsions d'un 
corps cyclique. Avant de laisser ce sujet nous 
voulons indiquer en peu de mots comment 
on peut traiter cette question dans le cas général d'un corps cyclique formé par 
un nombre quelconque de verges droites ou courbes très-minces soudées à leurs 
bouts en plusieurs nœuds. 

Soit AB l'une des verges, les efforts dans toutes les 
sections de la verge seront égaux. Prenons l'origine comme 
centre de réduction, et désignons par 

(ab) Y (ab) y{ab) v (ab) v (ab) Y (ab) 





vyu-") VK"-") VK"-") V\ uu ) V\ uu ) V 
vYj Ar, vV, y\ A vV K j\, 
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les caractéristiques des efforts correspondantes aux sec- 
tions de la verge. 



Soient x\ 



w J a ) J a ) J°) J a ) J a ) 



ry v ' ry *> ' ry \ • ry 



les composantes de la 



translation et de la rotation qui ont conduit la molécule A 
de l'état naturel à l'état actuel et x\ 



(b) Jb) Jb) Jb) Jb) Jb) 



ry v ' ry v ' ry v 
2 S 4 



% % 



les quantités analogues relatives à la molécule B. 

Soient e\ e 2 e£ e\ e 5 e£ les caractéristiques de la distorsion qu'on a 
appliquée à une section de la verge AB. On aura des relations linéaires 



.(a) 



x 



(b) 



,(«&) 



/ j H r iXi 



1, 2, 3, 4, 5, 6 



et les coefficients constants H ri dépendront de la forme initiale et de la constitution 
de la verge. 
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En même temps si nous prenons toutes les verges qui aboutissent au point A 
et qui sont soudées ensemble dans ce point on doit avoir pour l'équilibre 

x (ab) + x (ac) + x (ad) + . . . = Q , = 1, 2 , 3, 4, 5, 6. 

On a donc six équations pour chaque nœud et six équations pour chaque verge. 

Ces équations sont tout à fait semblables aux équations de Kirchhoff pour 
les courants électriques dans les fils ; seulement elles sont sextuplées. On peut 
développer une théorie du même type que celle de Kirchhoff où les caractéristiques 
des distorsions remplacent les forces électromotrices et les caractéristiques des 
efforts les intensités des courants. 
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leçon. 



1. Nous avons étudié dans les leçons précédentes les solutions des équations 
de l'élasticité qui appartiennent au type elliptique, et nous avons envisagé parti- 
culièrement les solutions polydromes et le rôle qu'elles jouent dans la théorie. Le 
cas de l'élasticité nous offre le meilleur exemple pour trouver leur interprétation 
et leur utilité pratique. 

Nous voulons pénétrer aujourd'hui plus profondément dans l'étude des so- 
lutions des équations du type elliptique, mais pour cela nous laisserons de côté 
le cas de l'élasticité, et nous envisagerons le cas le plus simple possible, celui de 
l'équation de Laplace qui correspond à la théorie du potentiel. 

On peut considérer l'équation de Laplace relative à un nombre quelconque de 
variables 



d 2 u d 2 u d 2 u d 2 u d 2 u d 2 u d 2 u d 2 u d 2 u 

dx 2 dy 2 dx 2 dy 2 dz 2 dx\ dx\ dx\ dx\ 



Le cas devient de plus en plus compliqué, mais bien des théorèmes s'étendent 
sans aucune difficulté de proche en proche. Par exemple le théorème de Green, 
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avec toutes ses conséquences et toutes ses applications. Il est utile, je puis dire 
même indispensable, pour cette étude d'introduire le langage emprunté à la théorie 
des espaces à plusieurs dimensions. C'est seulement un langage, mais il est un 
instrument très-puissant pour énoncer des théorèmes, et il amène par l'induction 
à trouver des propositions nouvelles. Je citerai dans cette direction les études de 
Ch. Bjerknes qui correspondent au mouvement d'une ellipsoïde à n dimensions 
dans un fluide à n dimensions. 

2. Cependant dès les premiers pas dans les recherches dont je viens de parler 
on a l'impression qu'il y a une lacune à combler, c'est à dire que quelque chose 
manque. 

Je vais citer un seul exemple. Les domaines à deux dimensions ont une seule 
sorte de connexion ; les domaines à trois dimensions ont deux sortes de connexions. 
En effet nous avons déjà distingué les espaces cycliques et les espaces acycliques, 
et nous avons vu quel rôle joue la cyclicité dans le cas de l'élasticité. Mais prenons 
le volume compris entre deux sphères concentriques. Il est acyclique, c'est à dire 
que tout cycle formé par une ligne fermée peut se réduire aussi petit que l'on veut 
par une déformation continue, sans sortir de l'espace, mais la même propriété ne se 
vérifie pas pour toute surface fermée, car une surface menée entre les deux sphères 
ne peut pas être réduite aussi petite que l'on veut sans sortir de l'espace renfermé 
entre les deux sphères. Il y a donc deux sortes de connexions dans le cas de trois 
dimensions, on les appelle connexion linéaire ou cyclicité et connexion superfi- 
cielle, mais on connaît une seule sorte de polydromie. Si le nombre des dimensions 
augmente, on a plusieurs sortes de connexions toujours plus compliquées, mais à 
cela ne correspondent pas plusieurs sortes de polydromies. 

Le cas de deux variables est étroitement lié à la théorie des fonctions. En 
effet par la séparation de la partie imaginaire dans une fonction d'une variable 
complexe on trouve deux fonctions de deux variables qui satisfont à l'équation de 
Laplace. Réciproquement toute fonction de deux variables qui satisfait à l'équation 
de Laplace possède une fonction conjuguée, et en ajoutant la première à la fonction 
conjuguée multipliée par i = V^ï on trouve une fonction d'une variable complexe. 

L'existence de la fonction conjuguée et la liaison avec la théorie des fonctions 
engendrent un ensemble de propriétés dans le cas de deux variables qui manquent 
dans le cas d'un nombre plus grand de variables et ce sont ces propriétés pourtant 
qui conduisent dans le cas de deux variables aux propositions les plus importantes 
et les plus cachées de la théorie. 

On est amené par là à soupçonner que la lacune dont j'ai parlé tout à l'heure 
pourrait être comblée en étendant au cas de plusieurs variables la théorie des 
fonctions conjuguées. 

Je consacrerai cette leçon et la suivante à donner un aperçu de quelques 
recherches qu'on a faites à ce sujet qui conduisent à deux branches différentes 
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provenant cependant d'une même source. 

3. Rappelons les propriétés des fonctions conjuguées dans le cas de deux 
variables. 

Si la fonction u satisfait à l'équation de Laplace A 2 u = en posant 







— =X —=Y 
dx dy 


on trouve 


dX 

dy 


dY dX dY 
dx ' dx dy 


C'est pourquoi 







0. 



X dx + Y dy, Y dx — X dy 

sont des différentielles exactes. 

Envisageons les fonctions dans une aire a à connexion simple, où X, Y sont 
finies, continues et monodromes. Calculons 



u 



u 



Ai 



(X dx + Ydy), 



vo 



■Il 



+ / (Ydx-Xdy), 




les intégrales étant étendues à une ligne 5 qui en partant 
du point A aboutit au point Ai, où uq est la valeur de u 
dans le point A et i>o est une quantité constante. On a que 
les intégrales ne dépendent pas de la ligne 5 pourvu que A 
et Ai ne changent pas et u est la valeur de la fonction u 
dans le point Ai et v est la valeur de la fonction conjuguée 
dans le même point. L'aire étant à connexion simple si le 
point Ai coïncide avec A les intégrales sont nulles, c'est 
pourquoi « et « sont des fonctions monodromes, mais si 
l'aire a est à connexion multiple u et v peuvent être poly- 
dromes. 

La dérivée de u dans une direction quelconque £ est égale à la dérivée de v 
dans la direction normale rj, si les directions £, rj sont orientées l'une par rapport 
à l'autre comme x et y. 

En calculant u + iv = f, on obtient une fonction de la variable complexe 
x + iy. 

4. Passons maintenant au cas de trois variables. 

Soit u une fonction qui satisfait à l'équation de Laplace ; posons 



du 
dx 



X, 



du 
dy 



Y, 



du 
dz 
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on aura évidemment 








K ' dy dz 


dX dZ 
dz dx 


dY 
dx 




(2) 


dX dY dZ 
dx dy dz 


= 0. 





0; 



La propriété caractérisée par les équations (1) peut s'énoncer en disant que le 
curl du vecteur ayant pour composantes X, Y, Z est nul et la propriété caractérisée 
par l'équation (2) en disant que la divergence du vecteur est nulle. 

■5. Cela posé, soit S un domaine ayant une connexion simple linéaire et 
aussi une connexion superficielle simple. Toute ligne fermée sera le contour d'une 
surface renfermée dans le domaine et toute surface fermée sera le contour d'un 
volume aussi renfermé dans le domaine. C'est pourquoi en employant le théorème 
de Stokes on trouvera 

(3) !{X dx + Ydy + Zdz) = 

J s 

l'intégrale étant étendue à une ligne fermée 5, et en appliquant le théorème de 
Gauss on aura 

(4) / (X cos nx + Y cos ny + Z cos nz) do = 

où cette intégrale est étendue à une surface fermée quelconque. 

On tire de l'équation (3) que si l'on étend l'intégrale à une ligne s' qui en 
partant d'un point A aboutit à un point Ai la valeur de l'intégrale dépendra 
seulement des points A et Ai, et ne dépendra pas de la ligne comprise entre ces 
points. 

Si nous appelons u la valeur de u au point A on aura 



U = Uq 



(Xdx + Ydy + Z dz) 

J s' 



u étant la valeur de la fonction au point A\. Envisageons maintenant l'intégrale (4), 
mais supposons qu'elle ne soit pas étendue à une surface fermée, mais à une surface 
ouverte. Pour déterminer la direction de la normale n à la surface nous convenons 
que, le bord de la surface ayant une direction déterminée personnifiée par une 
poupée qui regarde la surface, n soit dirigée de droite à gauche de la poupée. Cela 
posé on tire tout de suite de l'équation (4) que la valeur de l'intégrale étendue à 
la surface ouverte ne dépend pas de la surface où l'on a calculé l'intégrale, mais 
de son bord. C'est pourquoi l'intégrale étendue à une ligne conduit à une fonction 
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des points de l'espace : c'est à dire à la fonction primitive, et de même l'intégrale 
étendue à une surface conduit à une fonction des lignes de l'espace. 

Il n'y a pas de difficulté à concevoir des fonctions de lignes comme on conçoit 
des fonctions ordinaires, c'est à dire les fonctions des points. En effet, si l'on regarde 
tous les points de l'espace, et si l'on suppose qu'à chaque point correspond la valeur 
d'une quantité, on a une fonction des points de l'espace, c'est à dire une fonction 
ordinaire de trois variables. Mais toutes les lignes forment aussi des éléments 
géométriques de l'espace, et l'on peut évidemment concevoir une quantité qui a 
une valeur correspondant à chaque ligne. Voilà donc une fonction des lignes. 

Celle que nous avons trouvée n'est 
qu'une fonction particulière. Elle correspond 
à des lignes fermées, elle est continue et elle a 
aussi une autre propriété intéressante. Soient 
s et s' deux lignes fermées qui ont un arc 
commun AB. Si on doit parcourir l'arc AB 
en directions contraires en supposant qu'il 
appartienne aux deux lignes, on pourra re- 
trancher AB et obtenir une courbe s" dont 
la direction même sera déterminée. On écrira 




s . 



Soit V la fonction des lignes que nous avons trouvée, et désignons par V, V, V" 
les valeurs correspondant aux lignes s, s', s" il est évident que 

V" = V + V. 



C'est pourquoi l'on peut appeler la fonction des lignes que nous avons trouvée 
une fonction de premier degré. 

6. Elle a été obtenue par un procédé d'intégration ; il est facile d'y appliquer 
un procédé inverse qu'on peut appeler un procédé de dérivation. 

Revenons à une fonction ordinaire des points de 
l'espace. Pour la dériver il suffit de déplacer le point 
indice de la fonction dans une certaine direction, de 
calculer le rapport entre la variation de la fonction et 
le déplacement, et d'en déterminer la limite lorsque le 
déplacement devient infiniment petit. 

De même prenons la fonction des lignes que nous 
avons trouvée. Déplaçons en V un arc l de la ligne, et 
cherchons la variation correspondante de la fonction. 
A cause des propriétés que nous avons trouvées cette 
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variation sera la valeur de la fonction correspondant à la ligne formée par les arcs 
V et l en supposant que l'on invertisse la direction dans laquelle on doit parcourir 
l'arc l. Elle pourra donc s'exprimer par 

W = (X cos nx + Y cos ny + Z cos nz) cE, 

E étant une surface ayant pour contour les lignes l et V ; n sa normale prise dans 
la direction qu'on a établie précédemment. 

Supposons maintenant que la surface S en décroissant infiniment d'une manière 
régulière tende vers un point M ; on aura 

W 

lim — - = X cos nx + Y cos ny + Z cos nz 

en prenant les valeurs de X, Y, Z qui correspondent au point M. 

W 
On peut appeler lim — - la dérivée de V par rapport à S, et l'on peut la 

dV s 

représenter par — — . Il est évident que le signe dont la dérivée est affectée n'est 

Û2-t 

déterminé que lorsqu'on connaît la direction positive de la normale à la surface E. 
Prenons la dérivée de V par rapport à une surface S normale à l'axe X en M, on 
aura 

De même, si Si et S 2 sont des surfaces normales aux axes Y et Z, on aura 

dEi ' dS 2 
c'est pourquoi l'on peut désigner X, Y, Z par les symboles 

dV dV dV 



d(yz) ' d(zx) ' d(xy) 

Revenons maintenant à la fonction primitive u. Nous avons 

du du du 

dx dy dz 



donc 



du dV du dV 


du 


dV 


dx d(yz) ' dy d(zx) ' 


dz 


d(xy) 


ou en général, n étant la normale à la surface S, 






du dV 






dn <9E 
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On tire de là qu'en dérivant la fonction des lignes V par rapport à une certaine 
surface et en dérivant la fonction u par rapport à la normale à cette surface on 
trouve la même valeur. Donc la fonction u et la fonction V ont une relation 
réciproque tout à fait analogue à celle des deux fonctions conjuguées dans le cas 
de deux variables. C'est pourquoi dans le cas de trois variables, pour arriver à une 
fonction conjuguée de la fonction potentielle, il faut introduire les fonctions des 
lignes. 

7. On pourrait penser que ce nouveau concept est tout à fait abstrait, mais 
nous voulons montrer bien au contraire qu'il correspond à des concepts réels et à 
des vues pratiques. 

Il suffit pour cela d'envisager un champ électromagnétique constant. Sup- 
posons que nous ayons un pôle magnétique unitaire qui se déplace en prenant 
toutes les positions possibles dans un certain domaine externe à toute masse mag- 
nétique, électrique et à tout courant électrique. Le potentiel du champ par rapport 
au pôle sera une fonction qui vérifie l'équation de Laplace, c'est à dire une fonction 
harmonique. 

Mais évidemment on pourra aussi penser avoir un courant électrique unitaire 
qui parcourt un circuit fermé et qui se déplace en prenant toutes les positions et 
les formes possibles dans le même domaine. Il est bien naturel que nous calculons 
le potentiel du champ par rapport au courant électrique. On aura par là une 
idée complète de la nature électromagnétique du champ, mais en même temps 
nous aurons envisagé une fonction des lignes fermées du domaine. Il est facile de 
démontrer par un calcul très simple que cette fonction des lignes, c'est à dire le 
potentiel du champ par rapport au courant, est justement la fonction conjuguée 
que nous avons introduite tout à l'heure. 

On a donc que le nouveau concept a une base réelle et se présente même 
comme un élément nécessaire à envisager dans les théories de la physique mathé- 
matique. 

8. Les fonctions conjuguées ordinaires qui se présentent dans le cas de deux 
variables peuvent s'obtenir en particularisant celles que nous venons de définir. 
Il suffit en effet d'envisager un champ électromagnétique à deux variables. Alors 
il faut remplacer le pôle magnétique par un pôle logarithmique et le courant 
électrique par un point tourbillon et le passage est déjà fait. 

Il y a aussi un autre cas particulier qui présente beaucoup d'intérêt. C'est 
celui des potentiels symétriques. Dans ce cas il suffit de considérer les valeurs 
de la fonction conjuguée correspondant aux lignes circulaires normales à l'axe 
de symétrie et ayant le centre sur l'axe même. On arrive par là à des fonctions 
ordinaires qu'on appelle les fonctions associées que Kirchhoff avait introduites en 
partant de concepts tout à fait différents. 
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9. Jusqu'à présent nous avons envisagé les fonctions des lignes qui s'obtien- 
nent comme des fonctions conjuguées des fonctions qui satisfont à l'équation de 
Laplace ; mais on peut généraliser la chose. Soient X, Y, Z trois fonctions quel- 
conques définies dans un espace S dont la connexion linéaire et la connexion 
superficielle sont simples, et supposons qu'elles soient finies, monodromes et con- 
tinues, et que leurs dérivées aussi soient finies, continues et monodromes. Si l'on 

dX dY dZ 
dx dy dz 

on pourra toujours par le procédé que nous avons in- 
diqué précédemment calculer une fonction des lignes de 
l'espace. C'est pourquoi la condition pour l'existence 
de la fonction des lignes est la condition (2). Les con- 
ditions (1) ne sont pas nécessaires pour son existence. 
Elles sont nécessaires pour qu'elle soit la fonction con- 
juguée d'une fonction harmonique. Or, l'espace ayant 
les deux connexions simples, on peut calculer la valeur 
de la fonction des lignes pour chaque ligne fermée de 
l'espace dont la direction est connue, et pour chaque 

ligne on trouve une seule valeur. Envisageons maintenant un espace S dont la 
connexion superficielle ne soit pas simple, par exemple l'espace compris entre 
deux sphères concentriques ai et a 2 . Conduisons une surface a 3 comprise entre les 
deux, et traçons une ligne fermée s qui partage <r 3 en deux parties a 3 et <7 3 '. 

Ayant fixé la direction de la ligne s on pourra calculer la valeur de la fonction 
correspondante à la ligne s de deux manières en formant 




(4) 



(X cos n'x + Y cos n'y + Z cos n! z) da' 3 



et 



(5) 



/ (X cos n"x + Y cos n"y + Z cos n"z) da 3 , 



n' étant la normale à <r 3 et n" la normale à a 3 . Les deux normales ont évidemment 
la direction opposée, c'est à dire que si l'une est dirigée vers l'intérieur de a 3 
l'autre est dirigée vers la partie externe. Faisons maintenant la différence des deux 
intégrales. On aura l'intégrale étendue à toute la surface <r 3 , et on devra prendre 
la normale dirigée constamment dans la même direction. Mais a 3 à elle seule ne 
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limite pas un volume renfermé dans S, c'est pourquoi la dernière intégrale pourra 
n'être pas nulle et par suite les intégrales (4) et (5) pourront être différentes. 

On voit donc que dans ce cas la valeur de la fonction des lignes n'est pas 
déterminée pour chaque ligne. 

On peut se faire une idée de la polydromie d'une fonction des lignes aussi 
d'une autre manière. 

Pour avoir la différence des valeurs de la fonction cor- 
respondante à deux lignes s, s', on peut déplacer la ligne s 
d'une manière continue de sorte qu'elle vienne coïncider 
en position et direction avec la ligne s'. Soit a la surface engendrée par ce dé- 
placement, V la valeur de la fonction correspondante à la ligne s' et V la valeur 
correspondante à la ligne s. On aura 





V'-V 



(X cos nx + Y cos ny + Z cos nz) da, 




n étant la normale à la surface a dirigée par rapport à la 
direction de s' de la manière que nous avons indiquée dans 
le § 5. Déplaçons maintenant s' en lui faisant engendrer un 
tube de manière que s' se rapproche indéfiniment de s jusqu'à 
coïncider avec s', a deviendra une surface fermée. Si le do- 
maine où nous envisageons les choses a la connexion super- 
ficielle simple, a sera le contour d'un espace interne au do- 
maine et par suite en appliquant le théorème de Gauss on 
aura que l'intégrale étendue à a sera nulle, d'où l'on tire que 
les deux valeurs V et V seront égales. Mais, si le domaine 
a une connexion superficielle multiple, a n'est pas toujours 
le contour d'un espace interne au domaine et par suite les 
deux valeurs peuvent résulter différentes. Dans ce cas donc, 
si une ligne fermée en décrivant un tube revient à sa position 
initiale, en prenant des valeurs de la fonction qui se suivent 
avec continuité, on peut retourner avec une valeur différente 
de celle du départ. Par conséquent la polydromie pour les 
fonctions des lignes ne tient pas à la cyclicité, c'est à dire à 
la connexion linéaire, mais à la connexion superficielle. 
10. Nous avons vu que, pour rendre acyclique un espace cyclique, il faut faire 
une coupure ou plusieurs coupures. Qu'est-ce qu'il faut faire pour transformer un 
espace qui a une connexion superficielle multiple en un espace qui a une connexion 
simple ? 
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Envisageons l'espace compris entre deux sphères 
concentriques. Imaginons un tube infiniment mince qui 
réunit la grande sphère à la sphère interne. Par ce tube 
la connexion superficielle devient simple. 

Il y a d'autres cas où il faut un plus grand nom- 
bre de tubes pour obtenir la réduction. L'ordre de la 
connexion superficielle est mesuré par le nombre des 
tubes plus 1. A chaque tube correspond une constante 
spéciale pour la fonction des lignes obtenue par inté- 
gration, comme à chaque coupure correspond une constante pour une fonction 
ordinaire obtenue par intégration dans un espace cyclique. 

Si l'on cherche la fonction conjuguée de la fonction harmonique \ log(x 2 + 

y 2 ) dans le cas de deux variables on trouve la fonction polydrome arc tg-. De 

x . 
même si l'on cherche la fonction des lignes conjuguée de la fonction harmonique 

-. on trouve une fonction des lignes polydrome. 

y {x 1 + y 2 + z 2 ) 

11. On a calculé une fonction en regardant chaque ligne comme le contour 
d'une surface. Si l'espace est cyclique il y a des lignes qui ne forment le contour 
d'aucune surface. On peut trouver les différences des valeurs correspondantes à 
ces lignes, mais nous n'insisterons pas ci-dessus. 

Il y a aussi une autre manière pour calculer des fonctions des lignes de premier 
degré. Nous dirons en peu de mots en quoi elle consiste. 

L'équation (2) étant satisfaite on peut trouver trois fonctions Z 1 ,Z 2 -,Z 3 de 
x,y,z telles que 

X 
Y 
Z 



0Z 2 


dZ 3 


dz 


dy ' 


dZ 3 


dZ ± 


dx 


dz 


dZ x 


dZ 2 



dy dx 



Alors, si s est une ligne qui forme le contour de la surface o ayant pour 
normale n, on a par le théorème de Stokes 

/ [X cos nx + Y cos ny + Z cos nz) ds = [Z\ dx + Z 2 dy + Z% dz) 

J a J s 

en choisissant convenablement la direction de la ligne 5. On voit par là que la 
fonction d'une ligne peut se calculer par une intégrale étendue à la ligne même. 
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Nous avons fait usage précédemment du langage des champs vectoriels ; on 
peut y revenir pour donner une image visible et frappante des concepts que nous 
avons introduits. 

Lorsqu'un champ scalaire est divisé en lames infiniment minces telles que 
chacune correspond au même accroissement de la quantité scalaire, la fonction 
scalaire dans un point quelconque sert à dénombrer les lames comprises entre ce 
point et un point fixe. 

Envisageons maintenant un champ solénoïdal divisé en tubes équivalents in- 
finiment minces. Dans ce cas, pour connaître la nature du champ, il faut savoir 
dénombrer les tubes qui passent à l'intérieur de toute ligne fermée. Voilà le con- 
cept de fonction de ligne de premier degré qui surgit. Si le champ est en même 
temps solénoïdal et scalaire, le nombre des lames comprises entre un point fixe et 
le point variable et le nombre des tubes compris dans une ligne variable forment 
des fonctions conjuguées. 

12. Nous nous sommes bornés jusqu'à présent au cas de trois variables ; mais 
on peut généraliser les concepts que nous venons d'exposer au cas d'un nombre 
quelconque de variables. 

Dans un espace à n dimensions on peut imaginer des points, des espaces à 
une dimension, des espaces à deux dimensions, des espaces à trois dimensions etc. 
et des espaces à n — 1 dimensions, c'est pourquoi l'on peut imaginer des fonctions 
des points, c'est à dire des fonctions ordinaires de n variables, des fonctions des 
espaces à une dimension, des fonctions des espaces à deux dimensions etc. et enfin 
des fonctions des espaces à n — 1 dimensions. 

Si ces espaces sont fermés et si nous définissons une direction des espaces, on 
peut concevoir les fonctions des espaces à r dimensions de premier degré, c'est à 
dire, si S r et S' r sont deux espaces fermés à r dimensions ayant une partie com- 
mune qui a deux directions différentes selon qu'on la regarde comme appartenant 
à ^ ou à S' r , retranchons cette partie et envisageons l'espace S" formé avec les 
parties résiduelles. La fonction qu'on envisage sera de premier degré si, les valeurs 
correspondantes à S r , S' r , S" étant V, V, V", on a 

V" = V' + V. 

Les fonctions de premier degré des espaces à r dimensions peuvent se calculer 
par des intégrations étendues à des espaces à r + 1 dimensions, et l'on peut former 
des dérivées de ces fonctions tout à fait analogues aux dérivées des fonctions des 
lignes que nous avons envisagées. 

De cette manière dans un espace à n dimensions on peut obtenir des fonctions 

i n — 1 . . . n 

conjuguées de — - — types différents, si n est un nombre impair et de — types, si 
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n est pair. En effet à une fonction des espaces à r dimensions (r < n — 1) on peut 
faire correspondre comme fonction conjuguée une fonction de n — r — 2 dimensions. 

Dans un espace à deux dimensions on a une seule sorte de fonctions con- 
juguées, c'est à dire que les fonctions conjuguées à des fonctions harmoniques des 
points sont aussi des fonctions des points. 

Dans un espace à 3 dimensions les fonctions conjuguées à des fonctions har- 
moniques des points sont des fonctions des lignes, comme nous venons de voir. 
Dans un espace à 4 dimensions les fonctions conjuguées à des fonctions har- 
moniques des points sont des fonctions des espaces à deux dimensions, mais l'on 
peut avoir aussi des fonctions conjuguées à des fonctions des espaces à une dimen- 
sion qui sont aussi des fonctions de la même nature. 

Dans un espace à 5 dimensions les fonctions harmoniques des points sont 
conjuguées à des fonctions des espaces à trois dimensions, et l'on peut avoir des 
fonctions des espaces à une dimension conjuguées à des fonctions des espaces à 
deux dimensions. Ainsi de suite avec la règle que nous avons exposée précédem- 
ment. 

On arrive par là à une théorie générale des fonctions conjuguées où les fonc- 
tions harmoniques ordinaires à n variables se présentent comme un cas très- 
particulier. Il est bon de remarquer que dans les hyperespaces d'un nombre n 

TL — 1 

pair de dimensions il y a des fonctions conjuguées d'ordre — - — qui correspon- 
dent aux mêmes hyperespaces. 

La connexion de plus en plus compliquée des espaces à n dimensions joue un 
rôle très important et est étroitement liée aux différentes espèces des polydromies 
des fonctions, comme nous avons pu constater dans le cas de trois variables. Ces 
polydromies se présentent comme une extension de celles qu'on a dans les cas des 
intégrales Abéliennes. 

13. On peut envisager ces recherches comme une généralisation de l'étude 
des fonctions conjuguées, c'est à dire de la partie réelle et de la partie imaginaire 
d'une fonction d'une variable complexe considérées séparément. 

Mais si nous voulons parvenir à une vraie extension de la théorie des fonctions, 
c'est à dire de l'ensemble des deux fonctions conjuguées, l'une étant multipliée 
par y/— 1, il faut changer le point de vue, mais les mêmes concepts que nous 
avons exposés en sont toujours la base, c'est pourquoi nous avons dit dans la 
leçon précédente, qu'il y avait deux branches d'études qui provenaient d'une même 
source. Comme j'ai annoncé dans la première leçon je n'abandonnerai pas ce sujet 
sans avoir donné aussi un aperçu de ces recherches qui sont liées à la théorie des 
fonctions d'un côté et aux équations qui paraissent dans la physique mathématique 
d'un autre côté. 
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•M' 14. Si nous envisageons un plan et deux quantités complexes 

F et $ qui sont des fonctions des points du plan, la condition 
pour que F soit une fonction de $ est qu'en déplaçant le point 

indice M dans une direction quelconque MM' et, en calculant le 

AF 
rapport — — des variations des deux fonctions, la limite du rapport, lorsque MM' 

devient infiniment petit, ne dépend pas de la direction du déplacement mais soit 
une nouvelle fonction du point indice. Cauchy employait dans ce cas le mot de 
fonction monogène. 

Or une extension n'est pas possible dans l'espace, car si nous posons la même 
condition que nous venons d'indiquer dans le cas où le point indice au lieu de se 
déplacer dans un domaine à deux dimensions peut se déplacer dans un domaine 
à trois dimensions, on trouve qu'elle conduit à la même relation entre F et $ et 
l'on trouve la même chose dans le cas d'un espace à plusieurs dimensions. 

Mais nous avons vu que dans les espaces à trois dimensions et à plusieurs di- 
mensions le concept des fonctions ordinaires, c'est à dire des fonctions des points, 
n'est qu'un cas particulier du concept général de fonctions, car il y a les fonc- 
tions des lignes et en général les fonctions des hyperespaces. Si nous appliquons 
maintenant la condition que nous avons posée tout à l'heure pour F et $, mais 
en supposant que les fonctions soient des fonctions du nouveau genre que nous 
avons introduites, toute difficulté disparaît, et l'on peut étendre à un espace d'un 
nombre quelconque de dimensions la liaison de monogénéité de Cauchy, c'est à 
dire on peut étendre la théorie des fonctions au cas d'un domaine quelconque en 
trouvant quelque chose de nouveau. Faute du concept des fonctions des lignes et 
des hyperespaces la base où fonder la théorie manquait complètement. D'autre 
part il fallait bien penser qu'il aurait été nécessaire de tomber là dessus, si l'on 
voulait étendre la théorie des fonctions, car si l'opération de l'intégration dans 
le cas des fonctions d'une variable ne fait pas sortir du domaine des fonctions 
de la même nature, l'intégration des fonctions de plusieurs variables, c'est à dire 
l'intégration multiple qui est étendue à des domaines à plusieurs dimensions dont 
les limites sont des lignes ou des surfaces etc., devait nous amener au nouveau 
concept que nous avons introduit. 

15. Envisageons donc dans un espace à trois dimensions deux fonctions com- 
plexes de premier degré des lignes, et soient F et $ leurs valeurs correspondantes à 
une ligne L. Déformons un arc AB de L, et désignons par AF et A$ les variations 
de F et de $. Si en diminuant indéfiniment la déformation et la distance entre B 
et le point fixe A, le rapport 

A$ 



AF 
tend vers une limite qui dépend seulement du point A, on pourra dire qu'entre 
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les variables complexes $ et F passe une liaison tout à fait analogue à celle de 
monogénéité de Cauchy. 

On peut exprimer cette condition d'une autre manière en faisant usage des 
notations que nous avons introduites. Calculons 

d$ dF 

(6) ïë : dL = L 

E étant une surface quelconque, / doit avoir une valeur indépendante de la direc- 
tion de la surface S, c'est à dire doit être une fonction des points de l'espace. 

Décomposons maintenant $ et F et leurs dérivées dans les parties réelles et 
imaginaires en posant 

F = F 1 + iF 2 , $ = $i+i$ 2 , i = \/ r ï; 

dF dF dF 

P1 + IP2, 37 r = 9i + *92, 37 r = »"i + ir 2 ; 



d(y,z) ' " d(z,x) '" d(x,y) 

d$ d$ d$ 

u>i + iu> 2 , 37 r = Xi + %, 37 r=Pi+ïp2! 



d(y,z) ' d(z,x) ' d(x,y) 

la condition (6) pourra s'écrire 

(7 , u\ + JV2 = X1 + ÉX2 = Pi + ^P2 = , 

V '' Pi + *P2 çi + «92 ri + ir 2 

et l'on tire de là que $1 et $2 doivent vérifier les mêmes équations, c'est à dire 

d^ d%> d^ 

d(y,z) d(z,x) d(x,y) 

d^ d$ d^ d%> 

, . d ( l2 d(x,y) 13 d(z,x) \ d ( 23 d(y,z) 21 d(x,y) \ 

(8) ° = ^ D x ) + dy{ D 2 ) 

d$ d® 

'^31 37 7 — ^32" 



d I d(z,x) d(y,z) \ 

dz\ d, y 

où ^ peut être remplacé par $! et $ 2 et l'on a 

£11 =PÏ+pI, £22 = qf + ql, £33 = r\ + rj, 

£23 = £32 = qin + g 2 ^2, £31 = £13 = npi + r 2 p 2 , £12 = £21 = Pi<?i + P2Ç2, 

-Di = 92^1 - r 2 çi, D 2 = r 2 pi - p 2 ri, D 3 = p 2 q x - q 2 p 1 . 
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L'équation (8) remplace dans ce cas l'équation de Laplace et a bien des pro- 
priétés communes avec elle. 

16. On peut écrire l'équation (6) sous la forme 



(9) 



dF 



./'• 



Les fonctions $ et F peuvent s'appeler isogènes. Il est évident que deux fonctions 
des lignes qui sont isogènes avec une troisième sont isogènes entre elles. De même 
une fonction des points y? peut s'appeler aussi isogène avec F si 



(10) 



dF 
d(y,z) 



dx 



dF 



d(z,x) 



dy 



dF 



d(x,y) 



dz 



0. 



Elle sera isogène avec toutes les fonctions des lignes isogènes avec F. On démontre 
aussi facilement que / est isogène avec F, et réciproquement que si la fonction des 
points / est isogène avec la fonction des lignes F, on peut trouver une fonction $ 
telle que l'égalité (9) soit vérifiée. 

Toutes les fonctions des points qui sont isogènes avec F peuvent s'appeler 
isogènes entre elles. Puisqu'elles doivent vérifier l'équation (10) on en tire que deux 
d'entre elles sont indépendantes et les autres en sont des fonctions. Soient ipi et cp 2 
des fonctions indépendantes. A cause de l'équation (10) on pourra toujours écrire 



fil) 



dF 
d(y,z) 



dF 
d(z,x) 



dF 
d(x,y) 





d<pi 


d<pi 


<^3 


dy 

dip 2 


dz 
d(p 2 




dy 


dz 




d<pi 


dipi 


V?3 


dz 
dip 2 


dx 
d(p 2 




dz 


dx 




d(fi 


dipi 


^3 


dx 
dip 2 


dy 

dip 2 



x, 



y, 



dx dy 



où 9? 3 est une fonction des points isogène aux autres, c'est à dire tp 3 = (ps(<pi, ^2). 
Supposons que le domaine à trois dimensions que nous envisageons ait une 
connexion superficielle simple et que dans ce domaine les fonctions que nous con- 
sidérons n'aient pas de singularités, alors pour avoir la valeur de la fonction F 
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correspondante à une certaine ligne il suffira de regarder cette ligne comme le 
contour d'une surface o et de former 

F = (X cos nx + Y cos ny + Z cos nz) du. 

•la 

Soient u, v des coordonnées curvilignes de la surface. On aura 





dy 


dz 


du cos nx = 


du 
dy 


du 

dz 




dv 


dv 




dz 


dx 


du cos ny = 


du 

dz 


du 
dx 




dv 


dv 




dx 


dy 


do cos nz = 


du 
dx 


du 
dy 




dv 


dv 



dudv; 



dudv: 



dudv. 



C'est pourquoi l'intégrale précédente pourra s'écrire 



(12) 



F 



Vs 



dipi 


dfi 


du 


dv 


dip 2 


dip 2 


du 


dv 



dudv 



ip 3 dtpi dip 2 



On voit par là que tout ensemble de fonctions isogènes des lignes et des 
points dans un espace à trois dimensions ressort des fonctions de deux variables 
complexes et de leurs intégrales, comme tout ensemble de fonctions monogènes 
dans un espace à deux dimensions ressort des fonctions d'une variable complexe 
et de leurs intégrales. 

Si la surface a dans l'intégrale (12) est une surface fermée, l'intégrale devient 
nulle. Cette propriété n'est autre chose que la théorie de Cauchy pour les fonctions 
de deux variables telle que M. Poincaré l'a donnée. 

Si l'espace n'a pas une connexion superficielle simple, ou si les fonctions qui 
paraissent dans les formules ont des singularités, l'intégrale peut n'être pas nulle. 
Il faut alors remplacer la surface a par un ensemble de surfaces qui forment le 
contour complet d'un domaine à trois dimensions où les fonctions n'ont pas de 
singularités. 
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17. ifi étant une fonction de ipi et ip 2 telle que 



dtpi 



<P3(<Pl<P2), 



les équations (11) s'écrivent 



et l'équation (12) deviendra 



dF 
d(y,z) 



dF 
d(z,x) 



dF 
d(x,y) 



F 



d(p 4 


d(p 4 


dy 


dz 


d(p 2 


dip 2 


dy 


dz 


<9v? 4 


d(p 4 


dz 


dx 


dip 2 


dtp 2 


dz 


dx 


d(p4 


dp 4 


dx 


dy 


d<f 2 


d(p 2 



dx dy 



dif4dip 2 , 



d'où l'on tire que F peut s'obtenir par ip 4 et ip 2 ou qu'elle est composée par ces 
deux fonctions qu'on peut appeler des diviseurs de F. On peut écrire 

F = (^4^2)- 

La décomposition de F en ses diviseurs peut se faire d'une infinité de manières. 
Les diviseurs sont isogènes avec F et deux fonctions isogènes des lignes ont tou- 
jours un diviseur commun. Les diviseurs étant des fonctions des points, on peut 
les appeler élémentaires. 



'îeme 



leçon. 



18. La théorie qu'on vient d'exposer se rapporte au cas des fonctions de deux 
variables, et elle a l'avantage qu'on peut l'exposer sans recourir aux espaces à 
plus de trois dimensions. Mais si l'on se borne à ces espaces, on n'a pas évidem- 
ment une théorie complète ; et l'on voit la nécessité d'étendre tous les concepts à 
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l'espace à n dimensions. Il n'est pas difficile de le faire. Mais par là la théorie se 
développe d'une manière inattendue. C'est à dire qu'on ne trouve pas seulement la 
théorie des fonctions de plusieurs variables complexes et de leurs intégrales, mais 
quelque chose aussi de beaucoup plus général. En effet toute fonction complexe de 
premier degré d'une ligne peut être attachée par la liaison d'isogénéité à un en- 
semble d'autres fonctions des lignes et à des fonctions de deux variables complexes 
indépendantes, et peut se décomposer en facteurs élémentaires. La même chose 
ne se présente pas, lorsqu'on envisage une fonction des hyperespaces de r < n 
dimensions dans un hyperespace de n > 3 dimensions. Pour ces fonctions aussi il 
y a, comme dans le cas des fonctions des lignes dans l'espace à trois dimensions, 
une décomposition en fonctions d'un ordre plus simple qui présente des caractères 
arithmétiques analogues à la décomposition d'un nombre en facteurs premiers. 
Mais dans le cas des hyperespaces cette décomposition ne peut pas être poussée 
toujours jusqu'à arriver aux fonctions des points, c'est à dire que tous les diviseurs 
premiers de la fonction peuvent être aussi des fonctions des hyperespaces. Or, pour 
qu'une fonction soit liable à une autre par la liaison d'isogénéité il est nécessaire 
et suffisant qu'elle ait un diviseur élémentaire, c'est à dire un diviseur qui soit une 
fonction des points. Il faut donc poser des conditions spéciales avant de pouvoir 
établir la liaison d'isogénéité, et même après avoir établi cette liaison on a une 
théorie bien plus générale que celle qui ressort des fonctions de plusieurs variables 
et de leurs intégrales. Ce cas se rapporte seulement à l'ensemble des fonctions des 
hyperespaces dont tous les diviseurs sont des diviseurs élémentaires. 

Lorsqu'on envisage des espaces ayant un nombre de dimensions toujours crois- 
sant, les différentes sortes de connexions possibles augmentent toujours. Listing, 
Bei, Poincaré les ont étudiées. Elles jouent un rôle fondamental dans les théories 
précédentes, car elles déterminent les différentes sortes de polydromies des fonc- 
tions des hyperespaces. Dans le cas des fonctions d'une seule variable complexe 
on n'a qu'une seule sorte de polydromie qui ressort de la cyclicité, mais dans 
les cas précédents leur nombre est toujours plus élevé. La théorie des intégrales 
Abéliennes ne présente que le premier cas de cette théorie générale. 

19. Les fonctions des lignes et des hyperespaces nous conduisent naturelle- 
ment à dire quelques mots sur d'autres études très voisines qui ont bien des rap- 
ports avec des questions de la physique mathématique. 

Soient x,y,z les coordonnées des points de l'espace. Une ligne est déterminée 
par deux équations y = f(x), z = <p(x). Donc une quantité qui dépend d'une ligne 
peut être envisagée comme une quantité qui dépend des fonctions f(x) et <p(x). 

Il est donc tout naturel d'envisager les quantités qui dépendent de toutes les 
valeurs qu'une fonction ou plusieurs fonctions prennent dans certains domaines. 
C'est un concept qui se présente dans un grand nombre de cas. Le plus simple est 
celui d'une intégrale définie. La valeur de l'intégrale définie entre certaines limites 
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dépend de toutes les valeurs de la fonction entre ces limites. La physique math- 
ématique nous offre les exemples les plus frappants. Par exemple la température 
en un point d'une lame qui est chauffée au bord dépend de toutes les valeurs de 
la température au bord de la lame. Les fonctions des lignes et des hyperespaces 
rentrent évidemment dans ce concept. 

Il est bien clair que lorsqu'on parle d'une quantité qui dépend de toutes les 
valeurs d'une ou de plusieurs variables on entend quelque chose de bien différent 
d'une fonction de fonction. 

20. La définition que nous venons de donner est analogue à la définition 
de fonction ordinaire due à Dirichlet, d'après laquelle on dit qu'une variable est 
fonction d'une autre variable si à chaque valeur de cette quantité, entre des limites 
données, correspond une valeur de la première quantité. Nous n'entrerons pas 
dans les discussions très-importantes, qui se rapportent à cette définition. Il est 
bien évident qu'on a été amené à la donner par le concept fondamental de loi 
physique. M. Boutroux tout récemment a approfondi la question d'une manière 
fort intéressante. Mais en posant certaines conditions on peut passer du concept 
de fonction tel que l'a posé Dirichlet à celui de fonction analytique. Il est inutile de 
rappeler ces conditions qui sont bien connues et qui se rapportent à la continuité, 
à l'existence des dérivées etc. 

On peut procéder de la même manière dans le cas des fonctions qui dépendent 
de toutes les valeurs d'une fonction ou de plusieurs fonctions. 

Envisageons le cas le plus simple, celui d'une quantité F qui dépend de toutes 
les valeurs d'une fonction continue f(x) définie pour toutes les valeurs de x com- 
prises entre a et b. Il n'y a pas de difficulté à étendre le concept de continuité à 
la quantité F. Supposons maintenant qu'on part d'une fonction initiale /(#), et 
qu'on la change en la remplaçant par f(x) + eip(x) où e est une quantité infiniment 
petite. On peut tâcher de calculer la variation de F . Sous certaines conditions la 
partie du premier ordre de cette variation peut s'exprimer par une intégrale définie 



/ ^(6)^(6)^6 

J a 



La fonction F'(Çi) joue le rôle de première dérivée. Elle est indépendante 
de <p(x), mais elle dépend en général de toutes les valeurs de f(x), c'est pourquoi 
on peut tâcher de trouver la variation de F'(£i), lorsqu'on remplace f(x) par 
f(x)+e(p(x). Si l'on néglige les parties qui sont infiniment petites d'ordre supérieur 
à e, sous certaines conditions, on trouve que cette variation est donnée par 



e I ^(6)^"(6,6)d6 

J a 
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:i3) 



(14) 



F" (^1,^2) joue le rôle de première dérivée de F'(£) et de seconde dérivée de F. 
Elle est indépendante de <p{x), mais dépend de toutes les valeurs de f(x). Elle est 
une fonction symétrique de £1 et £ 2 - 

On peut aussi procéder au calcul de la troisième dérivée qui s'exprime par 
une fonction F'" (Ci, Ç 2 , £3) symétrique et ainsi de suite. 

Cela posé, on peut se proposer de développer la valeur de F qui correspond à 
/(#) + £(p(x) dans une série de puissances de e. Sous certaines conditions qui sont 
tout à fait semblables à celles qu'on a pour la série ordinaire de Taylor on trouve 

b 1 rb rb 



F = Fn 



<p(ti)F'(ti)dti 



1 



F"(ti,t2MtiM&)dÇidÇ 2 



3! 

1 
n! 



F'"(Çi, 6, âMàMàMà) #1 #2 dÇ 3 



J a J a 



où FW(^i,^ • ■ -tin) est une fonction symétrique des n variables, £1,^2 •■•£n ) in- 
dépendante de <p(x). 

En prenant e = 1 on passe facilement à une série analogue à celle de Mac- 
Laurin. 

21. On pourrait prendre comme définition d'une quantité F qui dépend de la 
fonction (p(x) définie pour les valeurs x comprises entre a et b l'expression précé- 
dente où l'on fait e = 1, en supposant que la série soit convergente. Ce serait une 
définition analytique d'une fonction qui dépend d'une autre fonction, et l'on pour- 
rait en étudier le domaine de convergence. Les différents termes représenteraient 
des fonctions de différents degrés. 

Cependant de cette manière on n'envisagerait pas le cas le plus général. 

22. Lorsque F dépend de toutes les valeurs d'une fonction ou de plusieurs 
fonctions, que celles-ci sont inconnues, et qu'il faut les déterminer, étant données 
certaines propriétés de F, alors on tombe sur les problèmes les plus intéressants 
de la théorie. Le calcul des variations n'est que le premier exemple d'une recherche 
de ce genre. 

Mais prenons le cas qui correspond à la formule (13), et supposons que nous 
écrivions 

0tt(x) = e<p{x)a + £ [ <p(Ci)F'^i,x)dCi + \^ I ! F"(Ç 1 ,Ç 2 ,xMÇ l MH 2 )dÇ 1 dÇ 2 

Ja ^ Ja Ja 



+ ■■■+-.8' 



F in \Çi,& • ■ .£»,zM£i) • • • if(U) dCi • • • d&, 



a étant une quantité donnée et \&(x) étant définie pour les valeurs de x 
comprises entre a et b. 
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Le cas le plus simple est celui où F", F"', . . . F^ . . . sont nulles, c'est à dire 
où l'on a 

(15) W(x) = ecp(x)a + e f v?(6)^'(6, x) d£ i; 

J a 

c'est le cas de l'inversion des fonctions de premier degré qui est le cas envisagé 
par M. Fredholm et par moi même lorsque la limite supérieure est égale à x. 

Le problème d'invertir une fonction générale de la forme (14) par rapport à 
l'inversion de la formule (15) n'est que le problème général de l'inversion d'une 
fonction transcendante par rapport à l'inversion d'une fonction de premier degré. 
Nous voulons montrer que le cas général peut être résolu. 

23. En effet pour obtenir <p(x) exprimée par ^f(x), c'est à dire pour invertir la 
formule (14), il suffit de tâcher de développer £<p(x) dans une série de puissances 
de ■&. 

C'est pourquoi calculons 

(d(e<p)\ (d(ey>)\ / rf 2 M \ = (<P(e<p)\ (d 3 (ep)\ [d 3 (eip)\ 

V dû J# =0 \ dû ; ' V d& ) ê=Q \ 

On trouve 

( d[ey(x)] \ [ h (djeipfr)]- 

V dû ) n L \ de 



dd 2 



ty(x) 



d# 3 7, =0 V dû* J 
F'(Ç 1 ,x)dÇ 1 . 



Ja \ ^ u / 

Si nous pouvons invertir cette formule, ce qui arrive par exemple lorsque le 
déterminant n'est pas nul, on trouve 



(16) 



[ de[ip(x)\ 

dû 



ip(x) = $(x)f3 



En dérivant une fois encore on a 



#(?7i)$'(?7i,£)d?7i. 



d 2 [e(p(x)} 
ââ 2 



a 



dû 2 



F > (Ç 1 ,x)dÇ 1 



fiy^^{^^i{^^> *«. 



fMi)ll fd[e<pfa)] 
M J 



c'est à dire 



(17) 



d 2 [sip(x)] 



dû< 



a 
o 

b r-b 



rf 2 M6)] 



dû 2 



F'(6,x)d6 



po pu 

/ / F"(^,^,x) 

J a J a 



MMÉon [d[eip(i 2 )} 



a *» a 

b r b 



v 



d:Ô 



V 



M 



dÇi d& 



$i(?7i, V2, x)^> (771) ^(772) drji dr]2 



39 



en remplaçant SÉEïM) , SÉEïM) par les ïaleurs données par la for . 

{ diï ) Q { d-d ) 
mule (16). Si l'on invertit la dernière formule en employant une formule analogue 
à la formule (16) on trouve 



d 2 [e(f(x)} 
èd 2 



b ,-b 



J a J a 



$"(?7i, V2, a;)*(ryi)*(ry 2 ) d^ dr\ 2 . 



De même on trouve 

d 3 [e<p(x)} ^ '•'' '•'' '•'' 



&"(Vi, V2, r]3,x)^(rii)^{V2)^{V3) drji dî] 2 dt] 3 . 



J a J a J a 

et ainsi de suite, de sorte qu'on aura 



de 



d 



1 



d 2 M 



,f 



î 

3! 



/•s 1 

•& I ^(r] 1 )^'(r] l ,x)dr] l + - 



d 3 (e(p) 
dïï 3 

b r b 



# 3 H = £^(x) = &^>(x)(3 



r 



$"(?7i, V2, x)^f {rji) ^(772) *?i ^î?2 



rîf 



H! 



$ (n) (77! . . . r? n , x)#(?7i) • • • ^(rj n ) dq x --- drj n 



Il faudra établir les conditions de convergence de la série précédente et son 
domaine de validité. — On peut regarder $( n ) comme une fonction symétrique de 
771,772 ■ ■ .i] n . Nous n'entrerons pas dans les détails par rapport à la vérification de 
la formule que nos avons donnée ni sur la question de l'unicité de la solution. 

24. Je ne puis terminer cet argument sans toucher à la question des méth- 
odes qu'on appelle de Jacobi-Hamilton. Les mémoires de Hamilton publiés en 1834 
et 1835 ont été le point de départ d'un grand nombre de travaux qui sont parmi 
les plus beaux de la mécanique analytique. Jacobi est le premier qui en ait fait 
ressortir tout l'intérêt. Le principe d'abord limité aux questions mécaniques a été 
étendu aux problèmes des isopérimètres et enfin à tous les problèmes du calcul des 
variations relatifs aux intégrales simples. On peut regarder la plupart des prob- 
lèmes de la physique mathématique comme ressortant des questions du calcul des 
variations relatives aux intégrales multiples, et il est par suite intéressant d'étendre 
les mêmes méthodes à ces cas. Pour atteindre le but il suffit de remarquer que la 
méthode Jacobi-Hamilton est fondée sur le concept d'envisager l'intégrale simple 
dont on veut annuler la variation comme une fonction des limites de l'intégrale et 
des valeurs arbitraires des fonctions inconnues aux limites. Donc pour étendre la 
même méthode au cas des intégrales multiples il faut les regarder comme dépen- 
dantes des lignes ou des surfaces ou des hyperespaces qui limitent les intégrales 
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données et de toutes les valeurs des fonctions inconnues au contour. C'est pourquoi 
les concepts que nous avons introduits sont nécessaires pour étendre les méthodes 
Jacobi-Hamilton aux problèmes de la physique mathématique. 
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1897. 

gième l eçon- 

1. Il est bien souvent nécessaire de revenir sur la signification des mots dont 
on fait un usage continuel dans la science. C'est leur sort d'évoluer et par suite de 
se transformer et de se plier à tant d'exigences différentes que non seulement le sens 
primitif est quelquefois presque complètement évanoui, mais que la signification 
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actuelle est très-difficile à découvrir lorsqu'on pense à la chercher. C'est ce qui 
arrive par rapport à la théorie générale des ondes. Tout le monde en parle et 
tous parlent de son influence, en l'opposant à d'autres systèmes. Mais les limites 
en sont elles-marquées, et peut-on établir sous quelles conditions une classe de 
phénomènes s'explique par le système des ondes ? Nous n'entrerons pas dans la 
discussion de ces questions qui mériteraient une étude approfondie ; nous dirons 
seulement qu'une grande classe de phénomènes peuvent se grouper au point de vue 
analytique sous une même catégorie bien définie, car leurs équations différentielles 
du type hyperbolique ont des propriétés communes qui ressortent principalement 
de l'existence de caractéristiques réelles. 

On envisage par là une grande partie des phénomènes qu'on a l'habitude de 
regarder comme dépendant du système des ondes, mais on ne peut pas dire que 
sous cette classe d'équations différentielles on ait groupé tous les phénomènes qui 
se relient à la théorie des ondes. On peut citer, par exemple, les ondes des liquides 
incompressibles et pesants qui dépendent d'équations d'une autre nature. (Voir 
ll eme leçon). Toujours est-il que l'étude des équations différentielles dont nous 
venons de parler forme une théorie que l'on peut envisager à part et approfondir 
par des méthodes générales qu'on est en train de perfectionner peu à peu. Par 
rapport à ces équations il faut dire qu'elles ont marqué le plus grand triomphe de 
la physique mathématique, car c'est l'analogie entre les équations des vibrations 
des corps élastiques et celles des champs électrodynamiques qui ont conduit à la 
théorie électromagnétique de la lumière et aux expériences de Hertz. 

2. Nous envisageons dans cette leçon quelques équations qui rentrent dans 
cette catégorie, et nous examinons les plus simples. Elles peuvent servir à nous 
montrer d'une manière claire, sans entrer dans les détails, les propriétés principales 
et la portée des méthodes qu'on a tâché d'employer pour leur étude. Considérons 
les équations 



d 2 v d 2 v d 2 v d 2 v d 2 v d 2 v d 2 v d 2 v d 2 v 

0, — + — - — = 0, — + — + — -7^2=0.... 



dx 2 dy 2 dx 2 dy 2 dz 2 dx\ dx\ dx\ dx\ 



La première est celle qu'on appelle des cordes vibrantes, la seconde est celle 
des membranes vibrantes, la troisième des corps vibrants etc. et l'on peut les com- 
parer avec les équations de Laplace à deux, trois, quatre variables indépendantes 
que nous avons déjà envisagées dans les leçons précédentes. On peut passer des 
unes aux autres par des transformations très-simples. C'est ainsi qu'en remplaçant 
dans la seconde par exemple z par iz, elle devient l'équation de Laplace à trois 
variables. Bien des propriétés peuvent se déduire par cette simple transformation. 

3. Nous voulons d'abord comparer deux formules qui nous montreront tout 
de suite le caractère des équations (1) et celui de l'équation de Laplace. Le cas 
le plus frappant et en même temps le plus simple est relatif à trois variables, et 



42 



l'on peut même considérer le cas plus général où le second membre n'est pas nul. 
Soit u une fonction qui vérifie l'équation 



(I) 



(f 



a 



d 2 u d' 



u 



f{x,y,z) 




(2) 



u(xi,yi,Z!) 



dx 2 dy 2 dz 2 

et supposons qu'elle soit finie, continue et monodrome 
dans un domaine S limité par un contour a . Supposons 
que la fonction / soit aussi finie, continue et mon- 
odrome. La formule de Green donne la valeur de w, 
dans un point M, de coordonnées £1,2/1,21 interne au 
domaine, exprimée par les valeurs de u et de la dérivée 
de u par rapport à la normale au contour. Voici la for- 
mule 

1 f f d Z du 1\ , 1 /" 1 



An 



II 



dn dn r 



da 



An 



dS 



a/(£i 



ou r = a/ (£1 — x) 2 + (2/1 — y) 2 + (z\ — z) 2 , x, y, z étant les coordonnées du point 
indice des intégrations. On suppose que la normale n soit dirigée vers la partie 
externe du domaine S. Si le point de coordonnées £1,2/1, z\ au lieu d'être un point 
interne Mj est un point externe M e , la formule précédente devient 

1 f ( d ~r 9U1\ 1 f ,1 „ 

(3) o = - [u-L- — )da-— f-dS. 

An J a \ on onr J Att J s r 

On peut par une simple transformation écrire la formule (2) d'une autre manière. 
Intégrons et dérivons successivement le second membre par rapport à Z\. On aura 



d 

dz\ 



1 

Att 



u 



dtp du 
dn dn 



'r 



da 



1 

Ait 



Z\ — z 



(*i - Z Y 



(2') u(xuVi,Zi) 

où l'on a posé 

* ~ log V \7(^i - x ) 2 + (yi - y) 2 ' V ( Xi - a; ) 2 + (^ - ^) 2 

de même la formule (3) pourra s'écrire 
(30 - 9 ' ' 



fipdS 







9^1 [ 47T 



u— - t^9? ) da — — I fipdS). 



dn dn 



Att 



s 



(II) 



4. Indiquons maintenant la formule analogue à la formule (2') pour l'équation 

d 2 v d 2 v d 2 v 
dx 2 dy 2 dz 2 
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Menons un cône de révo- 
lution à 45° ayant pour axe 
une droite parallèle à l'axe z et 
pour sommet le point M { dont 
les coordonnées sont X\y\Z\, 
et appelons a' une surface 
simplement connexe découpée 
par le cône et ayant pour 
contour la ligne d'intersection 
avec le cône. Soit S' l'espace 
compris entre le cône et la sur- 
face, n étant la normale à la surface dirigée vers l'intérieur de l'espace S', soient 
cosnx, cosny, cos nz les cosinus des angles qu'elle forme avec les axes coordonnés. 

Nous appelons (avec M. d'Adhémar) conormale la direction n' telle que 




cosn x 



cosnx, cosny 



cosny, cosn z = cosnz. 



On aura alors 



(5) v(x 1 y 1 z 1 ) 



d 

dz\ 



i r ( fy dv ; , 

V— — : — ^~i L P 1 d<J 



2n 



dn' dn' 



1 

2^ 



cp'fdS' 



ou 



(6) 



'r 



lot 



Z\ — z 



(zi - z)' 



y/( Xl - x) 2 + ( yi - y y y ^ ~ x ) 2 + fa* - y) 2 



en supposant v et f finies, continues et monodromes, et x,y,z étant les coordon- 
nées du point indice des intégrations. 

5. Les deux formules (2') et (5) sont tout à fait correspondantes l'une à 
l'autre. C'est en les comparant et en remarquant leurs différences qu'on pourra 
se faire une idée des analogies et en même temps des propriétés différentes des 
équations (I) et (II). D'abord on remarquera que dans la formule (2') la surface a 
et l'espace S renferment le point où l'on détermine la fonction u, tandis que dans 
la formule (5) paraissent seulement les valeurs de v et de la dérivée conormale le 
long de la surface a' découpée par le cône et les valeurs de / à l'intérieur de S'. 
Dans la formule (2') les espaces où l'on calcule les intégrales sont toujours les 
mêmes quel que soit M$. Dans la formule (5) ils changent en déplaçant M 4 . Les 
cônes de révolution à 45° ayant l'axe parallèle à z jouent donc dans ce cas un rôle 
considérable. Ils sont en effet les cônes caractéristiques. 
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En général lorsqu'on envisage une équation linéaire à coefficients constants 
de 2 ieme ordre 

v-^ v-^ ou v-^ , du , 



\J iXj fKJ *L> e 



OX,; 



r=l s=l i=l 

les surfaces ayant pour équations dans l'hyperespace x±, x-i . . . x n 



2—ir2—is^" rs v^ r •^rlv^s 



X 



0\ 



0, 



où x® • • . x° s sont des quantités constantes, s'appellent les cônes caractéristiques. 
Or par rapport à l'équation (I) les cônes ont pour équation 

(x - x ) 2 + (y- y ) 2 + (z- z ) 2 = 0, 

c'est pourquoi ils sont imaginaires, tandis que pour l'équation (II) les cônes sont 
représentés par 



(x - x ) + (y - y ) ~(z- z y 



0. 



Ils sont réels et sont justement les cônes de révolution à 45° que nous avons 
envisagés. 

6. Une autre particularité qu'il faut mettre en évidence est que dans la 

du 
formule (2') les valeurs u et — — au contour ne sont pas indépendantes entre elles. 

on du 

Si les valeurs de u sont données, les valeurs de -— sont déterminées aussi. 

on 

La chose est bien différente dans le cas de la formule (5). Pour certaines sur- 

dv 

faces a les valeurs de v et de — — sont indépendantes, tandis que pour d'autres sur- 
an 

faces ou certaines parties d'une surface elles sont indépendantes et pour d'autres 

ne le sont pas. Nous aurons l'occasion de revenir sur ce point. 

7. On peut chercher une formule qui 
corresponde à la formule (3') et il est bien 
facile de la trouver. 

Soit a' une surface interne au cône, et 
soit S' un domaine compris entre a' et la 
surface du cône, mais supposons que le som- 
met M e du cône soit un point externe au 
domaine S'. 

On aura alors 




(7) 



d 
dz\ 



1 

2^ 



», 



dip' du 
dn' dn 



~,V 



do 



1 

2^ 



ip'f dS' 



S' 
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(5') 



où n' représente toujours la conormale de la 
surface a'. 

8. En calculant la dérivée par rapport 
à Z\ qui paraît dans le second membre de 
l'équation (2') on revient à la formule origi- 
naire (2). 

Le calcul analogue dans la formule (5) 
n'est pas aussi simple, parce qu'une varia- 
tion des coordonnées x\,y\,z\ du point Mi 
change la surface a' découpée par le cône 

caractéristique. C'est pourquoi il est évident que la dérivation ne peut pas être 
faite directement sous les intégrales. La dérivée peut cependant se calculer très 
facilement dans le second et dans le troisième terme du second membre et par là 
la formule (5) peut s'écrire de la manière suivante 




v(xi,yi,zi) 



+ 



J__d_ 

2tt dzi } a , y/fa - zf - fa - x) 2 - fa - y) 2 

ci îs n : — — : -r [(x — x\) cos nx + {y — y\) cos ny\ ) v da' 

dv 



fa 



+ (2/i 
1 



y) 2 



da 



(yi 



y) 2 dri 

fdS'. 



1 

2~7T J ff , y/fa - zf - fa 
1 f 1 

2tt J s , y/ fa - z ) 2 - fa - x) 2 - (yi - y) 

La formule (7) peut se transformer de même. Dans ce cas le premier membre 
est nul. 

9. Revenons maintenant à la formule (2). Le second membre peut se dé- 
composer en deux termes. La première intégrale étendue à la surface a exprime 
évidemment une fonction harmonique, c'est à dire qui vérifie l'équation de Laplace 
à l'intérieur de S, c'est pourquoi si l'on prend 



on a 



WfayxZi) 



d 2 W d 2 W d 2 W 



dx\ 



fdS 



-4-îrf(x 1 y 1 z 1 ), 



dy\ dz\ 

formule qui correspond à un théorème célèbre de Poisson. De même prenons la 
formule (5') et considérons la dernière partie du second membre, c'est à dire 



W'fay x z x ) 



I. 



fdS 



s a/O^i - x ) 2 - (yi - y) 2 - ( z i - z ) 2 
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Puisque, comme on le vérifie facilement, la première partie du second membre 
satisfait à l'équation (II) où l'on suppose nul le second membre, on voit tout de 
suite que 



d 2 W d 2 W d 2 W 



dx{ dy( 



dz( 



-27rf(xiy 1 z 1 ) 



On a par là un théorème qui a bien des rapports avec celui de Poisson. Dans 
le cas de Poisson l'intégrale doit être étendue toujours au même espace S, dans le 
théorème que nous venons d'énoncer l'intégrale est étendue à la partie d'un espace 
qui est intérieur au cône caractéristique ayant pour sommet le point Xiyi z\. 

Par la formule (2) on peut facilement déduire les discontinuités relatives aux 
dérivées des potentiels des surfaces et les discontinuités des potentiels des doubles 
couches. On obtient des formules semblables en employant la formule (5') d'une 
manière analogue. 

10. Nous avons jusqu'à présent envisagé l'espace interne au cône, mais pour 
l'espace externe au cône on peut établir de nouvelles formules qui n'ont pas de 
correspondantes dans le cas de l'équation de Laplace. 

Menons le cône caractéristique ayant pour sommet le point M de coordonnées 
XiyiZi. 

Soit a" une surface externe au cône 
et découpée par le cône, et supposons que 
S" soit un domaine externe au cône limité 
par a" et par la surface du cône, n étant 
la normale à la surface a" dirigée vers l'in- 
térieur de S", appelons n" la conormalc. 
ir N n x ( , i ir Posons pour simplifier 




p = V( x ~ x i) 2 + (y-yi) 2 
v-vi 



cos np = cos nx 

P 



P 



cos ny. 



Enfin supposons que dans le domaine S", v et f soient finies, continues et 
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monodromes. Alors on aura 

1 



(8) 



v{x 1 y 1 z 1 ) 



cos np 



v do" 



27T 2 J a „ pv /( X - Xl )2 + (y- yi )2 _ ( 2 _ Zl )2 
1 d f 1 

Wdz~i h» y/{x-x x y + {y-y x y-{z-l^ 

((x-xxY + iy-yi) 2 - {z-zxY 
•log 



cos nz cos np ) v do 



2tt 2 J a „ y/( x - Xl )2 + (jj _ yi )2 _( z _ Zl y 

lo / (x - Xl ) 2 + (y - y t ) 2 - (z - z ± ) 2 \ dv ^„ 



P 



dn h 



2vr 2 J s „ ^(x - Xl )2 + ( y _ yi )2 _ ( 2 _ 2l )2 
' (x - Xi) 2 + (y - ?/i) 2 - (2 - ^) 2 



log 



fdS", 



où x,y,z sont les coordonnées du point indice des intégrations. Si nous écrivons 
dans le cas de l'espace externe une expression tout à fait analogue à celle qui 
paraît dans le second membre de l'équation (5) ou (5'), on ne trouve pas la valeur 
de v au sommet du cône, mais une valeur nulle. On a donc : 











1 



(9) 



dzi J a „ ^/(x - x x ) 2 + (y- y^ 2 - (z - z x ) 2 

[(x — X\) cos nx + (y — t/i) cos ny] ) v do" 

dv 



cos nz 



(xi - x) 2 + (y 1 - y) 2 
1 



do" 



a „ ^ _ Xl ) 2 + (y- y x ) 2 -(z- Zl ) 2 On" 

- 1 ^=f dS". 

s» a/(x - xi) 2 + (y - t/1) 2 - (z - 21) 2 

11. Particularisons maintenant les formules générales. On peut commencer 
par supposer / = 0. Il est évident alors que dans les formules précédentes tous 
les termes où paraît cette quantité s'évanouissent. Elles se rapportent alors aux 
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intégrales de l'équation 

(III) 



d 2 v d 2 



d 2 v 



dx 2 dy 2 dz* 



0. 



C'est pourquoi les formules (5') (9) (8) deviennent des formules correspondantes, 
mais d'une nature complètement différente, aux formules bien connues des fonc- 
tions harmoniques, c'est à dire des intégrales de l'équation de Laplace. Elles 
s'écriront alors 



vixiyiz-i) 



1 d 



(10) 



1 



27T d Zl J a , y /^ Zl - Z )2_ p 2 

1 /" 1 dv 



2vr J a , y /( 2l _ z )2 _ pi dn< 



cos nz cos np ) v do' 



do'. 



fil) 







dz! J a „ y/( Zl - z) 2 - p 2 
1 /" 1 dv 



cos nz 



Z\ — z 



cos np ) do" 



2vr J ' a „ ^( Zl - z )2 _ p 2 On' 



do" . 



(12) l 



v(xiyizi) 



cos np 



2tt 2 J a „ py J p 2 -[ z - Zl y 



-.v do" 



_1 d_ 

2^dz~i Ja » VP 2 " Tz - zi) 2 
1 f 1 



log 



p 2 - (z - Zlf 



cos nz 



Z\ — z 



cos np ) v do" 



27T 2 J a „ ^ p 2 _ ( z _ Zl )2 



log( p 2 -(.-^) 2 ^^ 



P 



<9n" 



12. On peut ensuite particulariser la forme des surfaces o' et o". Supposons 
d'abord que la surface o' soit le plan coordonné xy, et supposons z\ > z = 0. 
Alors la formule (10) devient 



(13) 



v(x 1 y 1 z 1 ) 



1 d 



1 



2^ 9*1 ^ v /( 2l _^)2_( :El _ ;r )2_ (2/l _ ?/) 2 

1 /" 1 dv 



v do' 



27T 7 ff , y /( Zl _ z )2_ {xi _ x) 2_ iyi _ y) 2 dz 

car la normale n et la conormale n' coïncident avec z. 



do', 
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Dans ce cas, v et — - sont indépendants entre eux 
oz 
(Voir § 6). Cette formule nous donne l'intégrale générale 

de l'équation (III) sous une forme bien connue qu'on 

appelle de Poisson-Parseval. 

On trouve une formule tout à fait analogue si l'on 
suppose que le point M* soit au-dessous du plan xy. 

Supposons maintenant que la surface a' devienne 
un cylindre indéfini 7 ayant les génératrices parallèles à 

l'axe z et aussi que la coordonnée z\ du point Mi soit supérieure aux coordonnées z 
des points de 7. Enfin supposons que v soit nulle pour toutes les valeurs de z 
inférieures à une certaine limite — Z. 

Dans ce cas la formule (10) deviendra 




(14) v{x 1 y\z 1 ) 




v(x,y,zi - 



<% 



Vë 



d_ 
dn 



v(x,y,zi -£) 



t% 



Vë 



où l'on suppose que 



P 



v^T 



r 



s étant l'intersection du cylindre 7 avec 

un plan parallèle au plan xy, et x, y étant 

les coordonnées des points de la ligne 5. 

Par rapport aux dérivées normales 

. , 1 1 d 5 

représentées par les symboles -— et — 

on on 
nous convenons que si / est une fonction 

de x,y,p et que si l'on regarde ces quan- 
tités comme des variables indépendantes 





s_l__ 

on 


df dx df dy 
dx dn dy dn 


Sf 


df dp 


on 


dp dn 



Dans ce cas les valeurs de v et de la 
dérivée normale de v sur le cylindre ne 
sont pas indépendantes entre elles. 

13. On peut enfin imaginer un cas 
mixte des deux précédents où la sur- 
face a' est constituée en partie par un 
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cylindre 7 et en partie par un plan parallèle au plan xy. Alors le second membre 
de la formule (10) est constitué de deux parties dont l'une a la forme du second 
membre de l'équation (14) et l'autre celle du second membre de l'équation (13). 
Sur 7 les valeurs de v et de la dérivée normale ne sont pas indépendantes, tandis 
qu'elles sont indépendantes sur le plan u>. 

14. Les formules (8) et (9) peuvent 
aussi se particulariser en supposant que a" 
se réduise à une surface cylindrique ayant les 
génératrices parallèles à l'axe z, et alors elles 
deviennent respectivement 



v(xi,yi,Z! -£) 





V7~ 



v(xi,yi,zi - 



<% 



^ 



(16) 



v(x 1 ,yi,zi) 




v(xi,yi,Z! -0 1o é 



S 
Sri 



v(xi,yi,zi -0 lo S 



P 2 -e 



p 2 -e 



dt 



-e 
dt 



Vp 2 ^ 



Vp 2 ^ 



dv 



Dans ce cas v et — - ne sont pas indépendants. 

on 



9 



îeme 



leçon. 



15. Dans le § 13 nous avons envisagé le cas mixte où sur u> les valeurs de v et 

dv 
de 7— sont indépendantes entre elles, tandis qu'elles ne sont pas indépendantes sur 

on 
la surface cylindrique 7. Sur cette surface il suffit de donner les valeurs de v ou de 

On dv 

—— pour déterminer v(x\ yi zA. Il faut donc tâcher d'éliminer les valeurs de — — ou 
on on 

de v sur 7. On peut atteindre le but dans certains cas très-aisément, et nous allons 

les approfondir pour montrer que la célèbre méthode des images découverte par 

Lord Kelvin pour l'intégration de l'équation de Laplace peut être aussi employée 

quelquefois pour les équations analogues ayant les cônes caractéristiques réels. 

On trouve par là un résultat complètement inattendu, c'est à dire que lorsqu'on 

envisage les équations à caractéristiques réelles, même dans le cas où l'on a un 

nombre infini d'images, la solution peut ne pas être donnée par des séries infinies, 

mais par un nombre fini de termes. Ce résultat semble au premier abord paradoxal, 
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mais il tient à ce que dans chaque cas le nombre des images dont on a besoin est 
fini en vertu de la réalité des cônes caractéristiques. 

Supposons que la surface cylin- 
drique 7 se réduise au plan yz (z > 0) 
et uj au plan xy (x > 0), c'est à dire 
que u> soit un segment d'un cercle et 
7 un segment d'hyperbole. 

Dans ce cas la formule (5') s'écrit 



(17) 



v{x 1 y 1 z 1 ) 




2ir dzi 



1 

2^ 



sf{z x - z) 2 - p 2 
1 dv 



cos(np)v d'y 



y/fa - Z) 2 - p 



2 dx 



d'y. 



dv 



Or pour éliminer les valeurs de v ou de — — sur 7 on peut se servir de la 
méthode des images. 




Soit M[ l'image optique du point M; réfléchi par le plan yz. Menons le cône 
caractéristique qui a pour sommet M[. Il découpera sur le plan yz le segment 
d'hyperbole 7 et sur la partie du plan xy (x > 0) un segment de cercle uj' . Si nous 
appliquons la formule (7), transformée convenablement, à l'ensemble des surfaces 
uj' et 7 on trouvera 

1 d [ 1 . , 1 f 1 dv 

2^ 



(18) 







2-7T ÔZ\ 

1 d 



Vzï 



du' 



P 1 



Vzï 



P 



2 dz 



du' 



1 



Zl 



2-7T dz\ 

-/— 



1 



z)'< 



- P 

dv 



ii 



cos(np)v d'y 



z) 2 - p 



2 dx 



d'y 
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ou 



a 



pi =(x 1 + xy + (y 1 - y y. 



En ajoutant les deux formules (17) et (18) on élimine v dans l'intégrale éten- 
due à 7, et par soustraction on élimine — — . La méthode des images conduit donc 

dx 
bien aisément au résultat qu'on cherchait. Passons maintenant au cas des images 

successives. Admettons que la surface 7 soit constituée par les segments d'hyper- 
boles 7' et 7" découpés par la nappe inférieure du cône ayant pour sommet Mi 
sur les demi-plans x = 0, x = a (a > 0) qui sont au-dessus du plan xy. Soit ou la 
partie de la bande du plan xy comprise entre les plans x = 0, x = a et renfermée 
à l'intérieur du cône. 

Dans ce cas la formule (2') s'écrit 



v{xiyizi) = 

1 d 
2tt dz\ 
_ 1 d 
2-k dz\ 



1 d 

2-k dz\ 



V4 



-v dui + 



VA 



p 



Ov 
2~d~z 



dio 



z\ 



VI 



Z\ 



1 



p 



cos{np)v d'y 



z\ 



VI 



z\ 



1 

2^ 

1 

cos(np)v d'y ' H 

p 2-k 



VI 



z\ 



Ov 



9 2 dx 



d'y' 



dv 



"VI 



z\ 



,2 dx 



d'y . 




x 



dv 



Dans cette formule les valeurs de v et de — sur 7 et 7 ne sont pas indépen- 

ox 
dantes, c'est pourquoi l'on peut tâcher d'employer la méthode des images pour 

éliminer les unes ou les autres. 

A cet effet prenons les images de Mi par rapport aux plans 7' et 7" et en- 
suite les images de ces images par rapport aux mêmes plans, et ainsi de suite 
indéfiniment. Le nombre des images qu'on obtient est évidemment infini, mais il 
est intéressant de remarquer qu'il n'est pas nécessaire de les employer toutes, car 
un nombre fini d'entre elles suffit pour la solution du problème. Menons les cônes 
caractéristiques dont les sommets sont les images que nous avons construites, et 
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envisageons les nappes inférieures de ces cônes, c'est à dire celles qui sont dirigées 
du côté négatif de l'axe z. Tous les cônes dont les sommets ont une distance des 
plans j' et 7" supérieure à z, ne rencontrent ces plans qu'au-dessous du plan u>, 

c'est pourquoi il n'est pas nécessaire de les envisager. Il est clair donc qu'on pourra 

dv 
éliminer les valeurs de v et de — sur 7' et 7" en faisant usage seulement des images 

ox 
qui sont à une distance de ces plans inférieure à z\. En employant donc dans ce 

cas la méthode des images, on ne trouvera pas des séries infinies, mais la solution 

du problème sera donnée par un nombre fini de termes. 

Au lieu des trois plans u, 7', 7" on peut en avoir cinq, dont quatre sont 

perpendiculaires au plan xy, et le rencontrent en formant les côtés d'un rectangle. 

La solution s'obtient toujours de la même manière par un nombre fini de termes. 

Si v s'annule sur les quatre plans perpendiculaires à xy, on obtient une solution du 

problème bien connu des vibrations d'une membrane plane, fixée dans un cadre 

rectangulaire sans avoir besoin de recourir à des séries, mais par des intégrales 

définies. 

16. On peut chercher à étendre les résultats que nous avons indiqués au cas 
des équations (1) à quatre, cinq etc. variables. Il suffit pour cela d'envisager des 
hyperespaces à 4, 5, . . . dimensions. Les cônes caractéristiques ont 3,4... dimen- 
sions. Au lieu des surfaces a', a" il faut envisager des espaces à 3, 4 . . . dimensions. 
Tout cela se conçoit très-bien. Mais si l'on passe de trois à quatre variables on 
arrive à un résultat inattendu, c'est à dire que les choses au lieu de se compliquer 
deviennent à un certain point de vue plus simples. En effet, on trouve des formules 
où la valeur de l'intégrale au sommet du cône dépend seulement des valeurs de la 
même fonction et de ses dérivées dans les intersections des cônes caractéristiques 
avec les espaces analogues aux surfaces a' et a". Il n'y a pas deux formules des 
types (2) et (8), mais une seule formule. 

En général, le cas de deux variables excepté, les choses se présentent d'une 
manière différente si l'on a un nombre pair ou si l'on a un nombre impair de 
variables, et le cas le plus simple est celui d'un nombre pair. Tout cela lorsque 
les équations ont une forme aussi simple que les équations (1). Si les équations 
avaient, par exemple, des termes de premier ordre la chose changerait. 

17. Nous n'entrerons pas dans les détails sur ce sujet qui sont très-compliqués. 
Nous renvoyons pour cela aux travaux de M. Tedone, de M. Coulon, de M. d'Ad- 
hémar et de M. Hadamard. Je rapporterai seulement les formules qu'on a dans le 
cas de quatre variables et qui ont été données pour la première fois par Kirchhoff, 
qui en a déduit avec rigueur le principe de Huyghens. 

Écrivons l'équation à quatre variables sous la forme 

d 2 v d 2 v d 2 v d 2 v 
v '' dt 2 dx 2 dy 2 dz 2 
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et envisageons x, y, z comme les coordonnées des points de l'espace à 3 dimensions. 
t peut être envisagée comme une quatrième variable, qui dans les questions de la 
physique mathématique est le temps. 

a étant le contour de l'espace S à l'intérieur duquel se trouve un point M\ 
de coordonnées £1,2/1, zi on a 

i 4-\ l f / d f v ( x ^y^ z ^- r )\ à fv(x 1 ,y 1 ,z 1 ,t-r)\\ 

4tt J a { on \ r J on \ r / J 

où r = y/(xi — x) 2 + (yi — y) 2 + {z\ — z) 2 , n est la normale interne à la surface a 

et les symboles — — et — - représentent par rapport à une fonction f(x,y,z,r) 
on on 

df df dr ôf df dx df dy df dz 
dn dr dn ' on dx dn dy dn dz dn 

en regardant x, y, z, r comme des variables indépendantes. Si la surface a se réduit 
à une sphère ayant pour centre le point M l5 la formule précédente se particularise 
et devient l'intégrale générale de l'équation (19) donnée par Poisson. 

18. Le succès de la méthode de Kirchhoff est dû à l'emploi des méthodes de 
Green et à l'existence de l'intégrale 

(20) ^M 

r 

de l'équation (19), où / est une fonction arbitraire. Pour toutes les équations pour 
lesquelles existent des intégrales de la même nature, qui n'ont d'autres points 
singuliers que r = 0, on peut procéder par des méthodes analogues à celles de 
Kirchhoff. 

Voyons si l'équation à trois variables qu'on peut écrire 

d 2 v d 2 v d 2 v 



dt 2 dx 2 dy 2 

possède une intégrale de la forme précédente, c'est à dire de la forme t?/(p±t) où 
û est une quantité indépendante de t, 

P= V( x i ~x) 2 + (2/1 -y) 2 
et / est une fonction arbitraire. En posant 

x\ — x = p cos a, y\ — y = p sin a 
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on trouve l'intégrale 

. , Asm \a + B cos \a , 

(21) 2 -— = 2 —f(p±t). 



Il paraît donc tout naturel que, même dans ce cas, on puisse employer cette 
intégrale en suivant les méthodes de Kirchhoff. Mais il est facile de se persuader que 
l'intégrale précédente est polydrome. Or les méthodes de Kirchhoff étant fondées 
sur celles de Green, c'est à dire sur le théorème de Gauss, dont nous avons parlé 
dans la deuxième leçon, on ne peut pas les employer dans ce cas. 

Il faut donc recourir à d'autres méthodes ou à d'autres intégrales dans le cas 
de l'équation à trois variables, et cela explique pourquoi les résultats étaient plus 
simples et aussi plus faciles à trouver dans le cas de quatre variables que dans 
celui de trois variables. C'est pourquoi le cas de 4 variables a été le premier à 
être résolu après celui tout à fait élémentaire de deux variables. Nous voyons quel 
rôle joue la polydromie des intégrales dans le cas des équations à caractéristiques 
réelles par l'exemple très simple que nous avons envisagé. 

19. Weierstrass a découvert une méthode pour donner l'intégrale générale 
des équations et des systèmes d'équations aux dérivées partielles à coefficients 
constants. Cette méthode a paru maintenant dans les œuvres de Weierstrass ; 
mais dès l'année 1881 il l'avait communiquée à M me Kowalewski qui chercha à 
l'employer pour intégrer les équations des vibrations de la lumière dans les milieux 
biréfringents. 

Supposons que 

(22) û(u,v,w) = t 

soit l'équation d'une surface fermée telle que toute droite qui part de l'origine ne 
la rencontre que dans un seul point, et supposons qu'on ait 

■&(Kx, Ky, Kz) = Kiï(x,y, z). 

Par 

f(u,v,w) du) 

représentons une intégrale étendue à l'espace renfermé entre deux surfaces (22) 
correspondantes aux valeurs t et t de t, dont la première est constante. 

Les formules fondamentales découvertes par Weierstrass sont les suivantes. 

Si 



(23) F(x, y, z, t) = D t \ ip(u, v,w)f(x + u,y + v, z + w) du 
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on a 



(24) 



de 



D X F = D t I <p—f(x + u,y + v,z + w)du-D t 



dip 
du 



f(x + u,y + v,z + w) du 



f d-â f dip 

D V F = Df I (p-—f(x + u,y + v,z + w) du — D t / -—f(x + u,y + v,z + w) du 
J dv J dv 

[d'à [dtp 

D Z F = D\ \ ip——f(x + u,y + v,z + w) du - D t j ^—f(x + u,y + v,z + w)du 



D f F = Dl I (pf(x + u,y + v, z + w) du 



où D a est le symbole de dérivation par rapport à a. On calcule de même les 
dérivées successives de F, et l'on tâche de prendre •& et y? de manière que les 
dérivées d'un certain ordre vérifient une relation linéaire à coefficients constants. 
On trouve alors une intégrale de l'équation ayant la forme (23) avec la fonction 
arbitraire /. De là on passe à l'intégrale générale. La méthode est analogue dans 
le cas des systèmes d'équations différentielles. 
Or envisageons la fonction 

ip(u, V, w, t) = cp(u, v, w)$(t — 1?) 

où $ est une fonction arbitraire. On aura les formules suivantes 

' dtp de . dip 



(24') 



du 

d%[) 

dv 

dip 

dw 

d'ip 

~dt 



dû 
-<p^-&(t-#) 



du 



$(*-#) 



du 



dv 
de 



dv 
dip 

dw 



$(£-??) 



Comparons les trois premières équations (24) et les trois premières équa- 
tions (24'). Les coefficients de / sous les intégrales des égalités (24) sont égaux et 
de signe contraire aux coefficients de $ et $' dans les égalités (24'). On ne trouve 
pas de changement de signe en comparant la dernière des équations (24) avec la 
dernière des équations (24'). 

On en déduit bien aisément la propriété suivante. Posons 



d n vb 



du hl dv h2 dw h z dt hi 



^ii hî 



d n F 



hnhi 



Q x h! Q y h 2 Q z h 3 gfh.4 



f h 1 h 2 h 3 h 4 
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où n = h\ + hi + h 3 + h^ on aura 

Fh 1 h 2 h 3 h i = D? +1 / P„/(x + u, y + v , z + w) du; + £>™ / P n -if(x + u,y + v, z + w) du 
+ ■■■ + D^ +1 I P h J(x + u,y + v,z + w) du. 

Tp hl h 2 h 3 h 4 = (-l) n - h HPn$ (n \t - 0) + Pn^-^it - 0) + • • • + P h M Hi) {t ~ ^)l 

et par suite si l'équation différentielle à coefficients constants 

( 25 ) Y. A hlh2h3h4 dxhi — — mh4 = 
possède une intégrale 

(26) F(x y zt) = D t / cp(u, v, w) f(x + u,y + v, z + w) dcj 

f étant une fonction arbitraire, l'équation 

Q n y 
( 25 ^(-l) n - fe4 A, lft2h3 , 4 — Qvh2 — mh4 = 

aura l'intégrale 

(26') ip(uv wt) = ip{uvw)<&{t-ti) 

où $ est une fonction arbitraire. Réciproquement si l'équation (25') possède l'in- 
tégrale (26'), $ étant une fonction arbitraire, l'équation (25) aura l'intégrale (26). 

L'intégrale (26') a la forme de l'intégrale (20), c'est pourquoi en comptant 
sur les singularités de ip pour u = v = w = 0, on peut passer de la méthode 
de Weierstrass à celle de Kirchhoff, c'est à dire des intégrales générales du type 
Poisson à des formules qui les comprennent. 

Le même résultat peut s'établir pour les systèmes d'équations différentielles. 

Je possédais en 1892 ce théorème lorsque j'étudiais le mémoire de M me Kowalewski 
sur la propagation de la lumière dans les cristaux. Puisqu'elle avait cru obtenir 
par la méthode de Weierstrass les intégrales générales des équations de l'optique, 
je croyais posséder la méthode pour parvenir aux formules analogues à celles de 
Kirchhoff et pour pouvoir établir le principe de Huyghens pour les milieux biréfrin- 
gents. Mais l'emploi de la méthode conduisait à des résultats paradoxaux. D'où 
provenaient-ils ? En voici la raison. Les fonctions qui dans ce cas remplacent <p dans 
les formules sont des fonctions polydromes du même genre que la fonction poly- 
drome que nous avons trouvée dans le § 18. Il n'était pas facile de s'en apercevoir, 
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car ces fonctions sont très compliquées, et elles avaient été étudiées depuis Lamé 
sans que cette propriété ressortît. C'est pourquoi ni la méthode de Weierstrass ni 
celle de Kirchhoff n'étaient applicables et les formules de M me Kowalewski ne don- 
naient pas les intégrales de l'optique. Voilà encore une fois que le rôle des fonctions 
polydromes ressort dans cette question d'analyse et de physique mathématique. 
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10 ième leçQn 

1. Nous avons jusqu'à présent examiné quelques équations du type ellip- 
tique et du type hyperbolique. Nous parvenons maintenant aux équations du type 
parabolique. 

L'équation de la propagation de la chaleur appartient à ce type. Fourier, 
Poisson, Laplace, pour ne citer que les plus anciens, ont consacré de profondes 
études à cette équation. Tout le monde sait que les méthodes introduites par 
Fourier sont devenues classiques et ont formé le point de départ, la base et le 
modèle de bien des recherches de l'analyse et de la physique mathématique. 

Plus récemment on a introduit dans la théorie de la propagation de la chaleur 
des méthodes qui se rattachent à celles de Green. Betti a consacré un très beau 
mémoire à cette étude, sans cependant recourir à la notion des caractéristiques. 
Appell et Boussinesq ont profondément examiné par des concepts plus modernes 
encore, bien des questions qui s'y rapportent. 

Mais à mon avis pour les équations du type parabolique on n'est pas si avancé 
que pour les équations des autres types que nous avons examinées dans les leçons 
précédentes. 

2. Nous allons envisager l'équation la plus simple possible 

du d 2 u 

Le rôle des deux variables t et z est très-différent, ce qu'on comprend dès le 
premier abord, car l'équation différentielle est de premier ordre par rapport à t et 
de deuxième ordre par rapport à z. 

Poisson a bien mis en lumière ce fait par un résultat sur lequel il a insisté 
beaucoup en y revenant plusieurs fois dans ses différents ouvrages. 

Nous allons l'exposer en empruntant à peu près les termes de son livre sur la 
théorie mathématique de la chaleur. 

Si nous essayons de développer u dans une série de puissances de t — h, h 
étant une quantité constante, nous trouvons, en calculant de proche en proche les 
coefficients, 

(2) u = V (z) + (t-h)-£ + -(t-h?^ + - + - l (t-hT 1 £ + - 
où ip représente une fonction arbitraire. Si l'on fait t = h, on a 

u = <p(z). 
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(p(z) est donc la fonction à laquelle se réduit u pour t = h. Au contraire si nous 
développons u par les puissances de z — k, on trouvera par la même méthode 

(z-k) 2 dip (z-k) 4 d 2 4> (z-k) 2n d n 4> 

I ~dt + 4! ~d¥ + '" + (2n)\ ~dF + '" 
' {z-kfdty (z-k) 2n+1 d n ^ 

3! ~dt + '" + (2n + l)! ~dF + " ' 

où -0(i) et \&(£) sont des fonctions arbitraires. 
En faisant z = k on trouve 

u = m, T^ = vl/ ( f ); 

du 
donc ?/>(£) et \&(£) représentent les valeurs de u et — pour z = k. 

Or il est facile de prouver que la même intégrale peut se mettre sous la 
forme (2) ou la forme (3). Il suffit pour cela de développer (p(z) dans une série de 
puissances de z — k et ip(t) et \I/(i) en des séries de puissances de t — h. 

En posant 

<p( z ) = a + ai (z -k) + ^-{z -k) 2 + --- + ^-(z-k) n + --- 

1-2 ni 

on trouvera en comparant les coefficients dans les formules (2) et (3) 

^(t) = a + a 2 (t -h) + ^-{t-h) 2 + --- + %i -h) n + --- 

II m 

$(t) = ai + a 3 (t -h) + ^(t-h) 2 + --- + °^(t - h) n + ■ ■ ■ . 

2! ni 

On voit donc que la même intégrale peut s'écrire avec une ou avec deux 
fonctions arbitraires de deux manières complètement équivalentes. 

Poisson a ensuite transformé les deux intégrales en les exprimant par des 
intégrales définies. Il écrit l'intégrale (2) sous la forme 



(2') 




u(z,t) 


1 

20F 


! 


v>(0 


e 4(t - /») 

Vt-h 


et l'intég 


;rale (3) sous la forme 










(3') 


u(z,t) - 


= 7 / i> ( 

J — oo 


v 4(c 


: + 




Vc + iÇ 






1 , 

+ 2^ - 


/•oo 

k) 

J —oo 


tf 


(- 


(z-k) 2 ^ 
4(c + iOv 
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ou 



1 

1 



•X 



dt 



. Vc + iÇ ' 
c étant une constante réelle qui n'est pas nulle. 

Après Poisson, Schlaefli est revenu là dessus en exprimant les mêmes intégrales 
sous des formes différentes dans un très beau mémoire inséré dans le T. 72 du 
Journal de Crelle. 

3. Cela posé voyons ce qu'on peut tirer des méthodes qui ont été les plus 
puissantes dans le cas des équations que nous avons envisagées dans les leçons 
précédentes. Il est évident que ces méthodes sont celles de Green, mais où l'on 
introduit la notion des caractéristiques par les méthodes de Riemann. 

Il n'y a pas de difficultés à procéder par une voie tout à fait analogue dans 
notre cas. 

On peut évidemment commencer par envisager une équation plus générale 
que l'équation (1), c'est à dire 



(!') 



d 2 



a 



dz' 



du 
~dt 



f(z,t). 



Les caractéristiques de l'équation (1') sont des lignes parallèles à l'axe z (voir 
la leçon 8 ième § 5). 

Prenons la direction de l'axe t de la gauche à la 
droite. Alors chaque caractéristique partage le plan 
en deux domaines qu'on peut appeler le domaine de 
gauche et celui de droite. 

Il est facile de démontrer le théorème : 
Soit s une courbe ouverte dont les deux bouts 
A, B se trouvent sur une caractéristique, et qui ren- 
ferme avec cette caractéristique une aire a située 
à gauche de celle-ci. Si nous connaissons u sur s 
et f(z,t) à l'intérieur de l'aire a , u sera déterminée 
* dans tous les points de a. 

La démonstration est très-simple et du même 
type des démonstrations des théorèmes analogues qu'on a dans les autres théories. 
En effet, si u\ et U2 sont des fonctions qui vérifient l'équation (1') et qui prennent 
les mêmes valeurs sur s, on calcule leur différence u% = u\ — U2 qui vérifiera 
l'équation (1) et qui sera nulle sur s. 

Par les transformations bien connues des intégrales on tire de là que 








u 3 



d 2 Us du3 

dz 2 dt 



da 



u\dl 



du 3 

dz 



da 
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où l est la partie de la caractéristique comprise entre les points A, B. 

C'est pourquoi u 3 sera nulle dans l'aire a et par suite U\ = u 2 - Donc il n'y a 
pas deux fonctions différentes qui satisfont aux conditions que nous avons posées, 
et on a par là la démonstration du théorème que nous avons énoncé. 

Il faut maintenant remarquer que les raisonnements que nous venons de faire 
ne pourraient pas être répétés si l'on supposait que la courbe s était à droite 
de la ligne caractéristique, et par suite il n'y a pas une symétrie des propriétés 
des intégrales de l'équation différentielle des deux côtés de chaque caractéristique, 
comme d'autre part il était facile de prévoir. 

4. On peut donner une formule tout à fait analogue à la formule de Green 
en envisageant avec l'équation (!') son adjointe 



(4) 



(5) 



d 2 v 

dz 2 



dv 
dt 



g(z,t). 



Alors on trouve 



(vf - ug) do 



du dv 



dz 



a 



dz 



cos nz 



uv cos nt 



ds 



uvdl, 



n étant la normale à la ligne s dirigée vers l'intérieur de o. 



Or la fonction qui joue le rôle principal dans 
cette théorie est, comme il est bien connu, l'in- 
tégrale suivante de l'équation (4) où l'on suppose 

{z-z{f 



2y/ïr^yt^t 



4(ti-t) 




t\ > t et z\ étant des valeurs constantes. Le radical 
est supposé positif. Prenons le point M 1 de coor- 
données ti Z\ à droite de la caractéristique l. Rem- 
plaçons cette valeur de v dans l'équation où l'on suppose g 
approcher M\ de la ligne caractéristique. 

A ce point il faut rappeler un théorème bien connu d'analyse. F(z) étant une 
fonction continue et finie et z" > z' on a 



0, puis faisons 



(I) 



lim 

t-L=t 



1 



y/h 



±{ti-t)F{z)dz 




Il faut prendre le second membre égal à 2yfnF{z-i) si z\ est comprise entre z' 
et z", et il faut le prendre égal si z\ n'est pas comprise entre les valeurs z' et z". 
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On suppose que le radical qui paraît sous l'intégrale soit positif. S'il était 
négatif le second membre changerait de signe. 

Deux cas pourront donc se présenter. Si le point M x s'approche d'un point M' 
compris entre A et B on aura 



lim 



tir 



dl 



égale à la valeur de la fonction u dans le point M'. Si au contraire il s'approche 
d'un point M" qui n'est pas compris entre A et -B la même limite est nulle. 
On tire de là les deux formules suivantes 




u cos nt > ds 



cos nz 



(6) 
(6') 

où t = ii est l'équation de la caractéristique l. Le premier cas se présente lorsque 
le point z 1 1\ est compris entre A et -B et le second cas lorsqu'il est externe. 

Dans cette formule, comme dans toutes les formules analogues que l'on trouve 
par les procédés de Green (Cf. le § 3 de la 8 lême leçon), la valeur de l'intégrale au 
point z\ t\ est donnée par les valeurs de l'intégrale même et de ses dérivées sur la 
ligne s. Or, comme nous avons vu (§ 3), il est nécessaire de connaître les seules 

valeurs de u sur l pour la déterminer dans a, c'est pourquoi il y a dans la formule 

du 
précédente la quantité — — qu'il faut éliminer. Elle s'élimine d'elle-même sur toute 

az 
partie CD de la ligne s qui est un segment de droite parallèle à l'axe z, car sur CD 

on a cos nz = 0. Ainsi, si les points A, B s'éloignent indéfiniment dans la direction 

positive et négative de z de manière que la ligne 5 devienne la parallèle t = h à 

ii i . , du , . . , 

1 axe z, la quantité — — est complètement éliminée. 

az 






s^- — 


A 


c 


f M' 






s 


l 


D 






s — 


B 
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Alors si nous posons dans la formule (6) / = elle devient la formule (2') de 
Poisson. En effet elle exprime la valeur de u par les valeurs de cette intégrale pour 
t = h. 



A, 



M"(z 2 t 1 ) 



t = h t = ti 




du 
5. Les valeurs de — sur s peuvent s'éliminer aussi en d'autres cas. Nous 
oz 
allons démontrer qu'elles peuvent s'éliminer sur toute partie CD de 5 qui est un 

segment de ligne droite telle que toute l'aire a est située d'un côté de cette droite. 

Prenons en effet le point M" symétrique du point M' par rapport au point E. 

M' étant interne au segment AB, M" sera externe. C'est pourquoi, si z\t\ sont 

les coordonnées de M 1 et z 2 h celles de M", en remplaçant dans la formule (6) z\ 

par Z2 on passe à la formule (6'). 

Le premier membre peut se décomposer en deux parties, l'intégrale qui est 

étendue à CD et toute la partie résiduelle. En appelant I cette dernière, on aura 



I' 



I" 



2nAF 



2nAF 



(zi-z) 2 
e 4(ti-t) 



CD 



y/t 



t 

(Z2 ~ Zf 

e 4(ti-t) 



CD 



Vh 



t 



du z\ — z 

u- 



dz 2(*i - t) 



du z 2 — z 

~dz ~ U 2{t 1 -t)_ 



cos nz — u cos nt \ ds = u(z\ ti) 



cos nz — u cos nt > ds = 



où /' et /" sont les valeurs de I si l'on se rapporte respectivement aux points M' 
et M". Ajoutons la première équation à la seconde après avoir multiplié celle-ci 
par une quantité constante — e c . Le second membre sera u{z\ t\) et dans le premier 
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du 
membre le coefficient de — — sous l'intégrale deviendra 

oz 



e 4 (*i 



t) 



tt+ c \ cosnz 
^ ; as. 



4(t 



(9-u 



Pour éliminer — — le long de 1 intégrale étendue à CD il suffit d annuler ce 
oz 
coefficient. Posons donc 

( Zl - zf (z 2 - zf 



4(ti - t) 4(ti - t) 
On tire de là en appelant £ e la coordonnée du point E 1 et h la longueur M' M" 

z e — z 2c 
*i - t ~~ h ' 

Mais la droite CD passe par le point E ; c'est pourquoi son équation est 



U -t 



a. 



a étant une quantité constante. 
Il suffit donc de prendre 



ha 



du 
dz 



-M' 



pour faire l'élimination de 

La méthode qu'on a suivie est bien une 
méthode des images, car on peut appeler le 
point M" l'image du point M' par rapport 
au point E. 

Si la ligne 5 est formée de plusieurs 
droites on peut appliquer dans certains cas 

le principe des images successives pour élim- 

du 
mer — — sur toute la ligne s. 
oz 
Soit s formée par trois droites AC, CD, 

DB dont CD est parallèle à z, alors la méth- 
ode peut servir pour tout point M' d'une manière très simple. L'on trouve dans 
ce cas un nombre infini d'images. 

6. Ce que nous avons exposé jusqu'ici nous prouve qu'une méthode tout à 
fait analogue à celle qu'on emploie dans le cas des équations du type hyperbolique 
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peut être aussi employée dans celui des équations du type parabolique, en tenant 
compte des lignes caractéristiques. 

Cependant nous pouvons bien facilement nous convaincre que par cette méth- 
ode on n'épuise pas les propriétés les plus importantes de l'équation (1'), et que 
la méthode même est impuissante à nous en dévoiler les plus cachées. 

En effet nous avons montré dès le 2 lème § le résultat obtenu par Poisson, 
c'est à dire que l'intégrale générale de l'équation (1) peut s'exprimer par une 
fonction arbitraire ou par deux fonctions arbitraires. Or la formule (6) peut nous 
conduire à l'intégrale (2), comme nous avons déjà vu, mais elle ne nous amène pas 
à l'intégrale (3). Même en tâchant de faire coïncider la ligne s avec l'axe des z on 
n'arrive pas à cette formule. 

On serait conduit à cette conclusion, que les équations du type parabolique 
n'offrent pas de prise aux méthodes modernes d'une manière complète, et que par 
exemple les anciennes méthodes de Poisson nous disent davantage. Il semble donc 
que les méthodes modernes peuvent nous donner l'explication d'un grand nombre 
de faits, mais tous les faits connus semblent ne pas ressortir d'elles. 

Malgré cela nous allons voir que cette impuissance n'est qu'apparente et que 
nous pouvons réunir dans un ensemble unique tous les faits et par une seule formule 
embrasser toutes les formules connues. Pour rejoindre le but il n'y a qu'à aborder 
les nouvelles méthodes avec une conception nouvelle, c'est à dire en les étendant 
aux variables complexes, ce qu'on n'avait pas fait jusqu'à présent. L'équation (1') 
se prête avec beaucoup de facilité à cette généralisation, car elle n'a que deux 
variables z et t, et si nous posons t = x + iy en regardant x,y,z comme les 
coordonnées des points de l'espace, on pourra se servir des images géométriques 
dans un domaine à trois dimensions. 

Envisageons une fois encore la fonction qui joue le rôle principal dans ces 
recherches, c'est à dire 



(zi - zf 



p. 4(ti-t) 

(7) 



Vû=t 



Nous verrons dans la leçon suivante que c'est la polydromie de cette fonction qui a 
l'importance la plus grande dans tous les résultats que nous allons obtenir. Cette 
fonction est en effet polydrome par rapport à la variable t. Or cette polydromie 
ne peut pas ressortir si on n'envisage pas t comme une variable complexe. C'est 
pourquoi toute la théorie que nous avons développée qui était bâtie en dehors de 
cette propriété était justement incomplète. 
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11 



îeme 



leçon. 



7. La première question qui se pose est de voir si l'on peut étendre dans les 
intégrales de l'équation (1') la variabilité de t, c'est à dire si l'on peut la regarder 
comme un nombre complexe. A cet effet considérons les intégrales de l'équation (1) 
qui est plus simple que l'autre. Il est facile de démontrer, et d'autre part la chose 
est connue, qu'elles ne sont pas en général des fonctions analytiques de t. 

Pour s'en persuader d'une manière tout à fait élémentaire prenons la fonction 



u(z,t) 



V(f)e 



4(t 



0(«-0" 5 # 



Si 



où </?(£) est une fonction finie et inté- 
grable. Si </?(£) est définie dans le do- 
maine (Oa) on voit aisément que u(z,t) 
est définie pour toutes les valeurs posi- 
tives de z à l'intérieur de la bande Si 
limitée par l'axe t, l'axe z et la droite 
t = a. 

On vérifie qu'à l'intérieur de ce do- 
maine u(z,t) est une fonction continue, 

, , . d 2 u du 

que les dérivées — — et — sont deter- 

oz 2 ot 
minées et que l'on a 



2a 



(1) 



d 2 u 
dz 2 



du 
~dt' 



Cela posé, si nous prolongeons </?(£) entre a et 2a, alors u(z,t) pourra se 
prolonger à l'intérieur de la bande S 2 comprise entre l'axe t et les droites t = a, 
t = 2a. 

Dans le domaine formé par l'ensemble des points internes aux bandes Si et 
E 2 u(z,t) sera continue et vérifiera l'équation (1). Mais puisque le prolongement 
de </?(£) peut se faire avec la seule condition qu'elle soit finie et intégrable et que 
d'autre part elle est tout à fait arbitraire, il s'en suit que les valeurs de u en E 2 
pourront être changées, celles en Si restant toujours les mêmes. 

Si par exemple nous prenons d'abord </?(£) = dans l'intervalle (a, 2a) et 
après </?(£) = 1, les valeurs de u(z,t) en S 2 seront changées. Cela prouve que 
u(z,t) n'est pas en général une fonction analytique de t. 

8. Il faudra donc poser comme une condition nouvelle que l'intégrale de 
l'équation (1) soit une fonction analytique de t, et de même si nous envisageons 
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l'équation (1'), il faudra poser deux conditions, c'est à dire que f et u soient des 
fonctions analytiques de t = x + iy. 

Supposons aussi, en ayant égard à l'équation adjointe (4), que g et v soient 
des fonctions analytiques de t = x + iy. 

Nous écrirons alors 



u(z,t) = u(x,y, z) 
v(z,t) =v(x,y,z) 



f{z,t) = f{x,y,z) 
9(z,t) = g(x,y,z) 



et l'on regardera x, y, z comme les coordonnées des points de l'espace. Il est évident 
que u, v, f,g seront en général des quantités complexes et on aura 



(8) 



du 
dt 



du 

dx 



.du 
' dy 



dv 
dt 



du 
dx 



.dv 
dy' 



X 



Posons 



du dv 

w â — M ^"' Y 
dz dz 



du dv 



u - 



dz dz 



Z 



ne 



et, en prenant une surface a ayant deux faces 
distinctes et pour contour la ligne s, appliquons 
le théorème de Stokes (2 ieme leçon § 2) en 
supposant que toutes les quantités précédentes 
soient finies, continues et monodromes sur la 
surface a. 
On aura 




G 




(9) 




fdZ _ dY\ fdX _ dZ\ 

\dy dz J \ dz dx J 



cos ny 



dY dX 
dx dy 



cos nz > do 



= (Xdx + Ydy + Zdz). 

J s 

Or on a bien aisément en vertu des relations (8) et des équations (1') et (4) 

—ivf + iug 
vf-ug 

0; 



dZ dY 

dy dz 

dX dZ 

dz dx 

dY dX 

dx dy 
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par suite l'égalité (9) deviendra 



(vf — ug)(cos nx + i cos ny) do 



du dv 

'— u— ) (dx + idy) + uvdz 

oz oz 



c'est à dire 



(10) 



\j} Vf - Ug) ta d(J 



du dv\ 
v- u— ) dt 

oz oz 



uvdz 



9. Cela posé supposons que dans tout le domaine à trois dimensions que nous 
envisageons u, f soient finies, continues et monodromes. Prenons 



fil) 



étant t 



4(t - t) 



Vk 



9 = 
t 




Si la droite x = xo, y = yo ren- 
contre la surface a, la formule (10) 
n'est pas applicable, car pour x = x , 
y = y la fonction v a une singu- 
larité. Mais elle n'est pas applicable 
en d'autres cas où la même droite 
ne rencontre pas la surface. En effet 
supposons que la surface soit limitée 
par deux lignes, et qu'en conséquence 
le contour 5 soit formé de deux par- 
ties, et supposons aussi que la droite 
x = xo, y = yo passe au milieu du 
bord interne de la surface, et que par 
suite on puisse parcourir sur la surface un cycle fermé c autour de la droite. Même 
dans ce cas la formule ne sera pas applicable, car v sera polydrome sur la surface 
à cause du radical qui paraît dans le dénominateur. En parcourant le cycle c et 
en prenant les valeurs de v qui se suivent avec continuité on reviendra au point 
de départ avec la valeur initiale changée de signe. 

Donc la formule (10) n'est applicable, en prenant v donné par la formule (11), 
que si aucun cycle fermé sur la surface n'entoure la droite x = Xo, y = yo et si en 
même temps cette droite ne rencontre pas la surface. 

Le contour étant formé par une seule ligne continue et fermée, il suffit évidem- 
ment de la seconde condition, pour que la première soit aussi satisfaite. 
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10. Prenons maintenant (fig. a) un con- 
tour s formé par une seule ligne continue, 
constituée par la courbe MNP et par la 
courbe M'N'P' reliées entre elles par les 
droites PP' et MM' parallèles à l'axe z. 

La ligne MNPP'N'M'M sera une ligne 
continue. Les flèches indiquent comment il 
faut la parcourir. Supposons que nous me- 
nions une surface o ayant pour contour cette 
ligne, et ensuite faisons approcher les deux 
bords MM' et PP' d'une même droite HH' 
du plan xz parallèle à l'axe z de manière 
qu'ils se disposent le long de cette droite. 
Mais en même temps supposons que si cer- 
taines parties des deux bords se superposent, 
comme il arrive dans la fig. b, les deux 
bords restent toujours distingués comme les 
bords d'une coupure, et qu'on doive les par- 
courir en sens opposé. La même chose est 
indiquée dans la fig. b' où les dispositions 
sont changées. Les flèches désignent toujours 
la manière de parcourir le contour. On peut 
donc regarder toujours la surface a comme 
limitée par un contour formé par la ligne 
continue et fermée MNPP'N'M'M. 

La ligne HH' ayant pour équation 



xi, y = 0, 



menons une droite parallèle KK' dans le 
plan xz ayant pour équation 



x 



x > xi, y = y = 



qui ne rencontre pas la surface a. Alors on 
pourra appliquer la formule (10), et l'on 
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trouvera en remarquant que g = 



MNP 



du ôva , 

va â w ^ — «* + UU A dz 

oz oz 



du dvA , , 

«An ô u ^r~ dt + WW A «Z 

OZ oz 



P'N'M' 



pp , 



1 



uv Ao dz + I uv Ao dz = - I v A J — da, 



M'M 



{z-zif 



e 4 ( x ° ~ *) 
où l'on a désigné par i^ la fonction — , A étant le point de coordonnées 

V x o - t 
x ,0,z 1 . 

11. On peut maintenant distinguer deux cas. Projetons la surface a sur le 
plan xy. Dans les figures nous avons désigné la projection par l'aire hachée. Le 
point xo, du plan xy peut se trouver entouré par la projection (fig. b) ou peut 
ne pas être entouré par elle (fig. b'). 

Puisque la fonction v est polydrome et qu'elle change de signe en parcourant 
un cycle autour de la droite x = x , y = 0, les choses seront différentes dans les 
deux cas. Les valeurs de v dans les points de la droite HH' sont réelles. Prenons 
ces valeurs positives sur P'P, alors 

dans le premier cas les valeurs de v seront négatives sur MM' 

dans le second cas les valeurs de v seront positives sur MM'. 

Or déplaçons le point A sur A Ai parallèle à l'axe x en l'approchant indéfin- 
iment du point Ai sur la droite HH', et supposons que la droite KK' se déplace 
avec A en restant parallèle à l'axe z, et supposons aussi qu'en se déplaçant elle 
ne rencontre jamais la surface a. 

Cherchons les limites de 



uva dz, 



PP , 



/ uv Aq dz. 

J M'M 



le 



Il suffit pour cela de faire usage de la formule (I) qu'on a déjà employée dans 
4. Puisque va est positif sur PP' on aura 



lim 



uva q dz = 2v / 7T'u(ti Z\), 



pp, 



si le point A x est compris entre les points P, P' et si le segment PP' a la direction 
positive de l'axe z. On aura au contraire 



lim / uvA Q dz = —2\/7ru(tiZi), 
! PP , 
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si le point Ai est compris entre les points P, P' et si le segment PP' a la direction 
négative de l'axe z. Enfin on aura 



lim / uva dz = 0, 



pp' 

si le point A\ est externe au segment PP' . (C'est le dernier cas qui se présente 
dans les figures b et b', mais on imagine tout de suite comment il faudrait dessiner 
les figures pour que les autres cas se présentent.) 

Les trois formules précédentes peuvent se réunir en une seule en représentant 
par le symbole (p) le numéro ±1 selon que le segment AiP a la direction positive 
ou négative de z. On aura alors 

lim / uv Ao dz = [(p') - (p)] y/ïr u(h z x ). 
J pp> 

Si nous étendons la signification du symbole (p), et si nous supposons qu'il 
représente 0, si Ai coïncide avec P, alors la formule précédente nous donne la 
limite de l'intégrale même dans le cas où Ai coïncide avec l'un des points P ou P'. 

Passons maintenant à 

lim / uvaq dz, 

J M'M 

si le point xq, est entouré par la projection (1 er cas) ; on a que va est négatif 
sur MM', c'est pourquoi 

lim / uva dz = — [(m) — (m')] y/ïr u(t\ z{). 

JM'M 

Si le point x , n'est pas entouré par la projection (2 lème cas) v Ao est positif 
sur MM' et par suite 

lim / uvA dz= [(m) — (m')~\\/ïru(ti Zi). 

JM'M 

Les deux formules peuvent se réunir en une seule en indiquant par le symbole 
(e) le numéro ±1 selon que la projection entoure ou n'entoure pas le point xo,0 
et l'on aura 

lim / uva dz = — (e) [(m) — (m)] \fïîu(ti z\). 

JM'M 

Donc 



lim 



uv Ao dz + I uva dz 

PP' JM'M 



{[(p')-(p)]+(e)[(m')-(m)]}^u(tiz 1 



d'où l'on tire en posant 
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[(PO " (P)] + (e) [(m') - (m)] 



/i 



(II) 



/iu(*i*i) = Uni ( -T / v A J — da - / 



MATP L 



P'N'M' 



du dv Ao 

du dv Aa 
va g u ~fj~ i dt + UVa ° dz 



Cette formule est la plus générale pos- 
sible et comprend toutes les autres. Il n'y a 
qu'à la particulariser pour déduire toutes les 
formules de la théorie. 

12. Supposons d'abord que le second 
cas se présente. Alors on peut très bien imag- 
iner que les lignes MNP et M'N'P' appar- 
tiennent au plan xz de manière que a soit 
une partie du plan, et l'on peut supposer 
que la ligne M'N'P' se rétrécisse indéfini- 
ment de sorte qu'elle se réduise à un point. 
Alors dans la formule (II) l'intégrale éten- 
due à M'N'P' s'annule, et elle conduit aux 
formules (6) et (6') dont la première nous 

amène à la formule (2') de Poisson. 

13. Soit M' N' P' une courbe appartenant au plan parallèle à xy mené par 

A , alors on démontre facilement que la limite de l'intégrale étendue à la ligne 

M' N' P' dans la formule (II) est nulle. En même temps on a (m') = (p') = par 

suite la formule (II) devient 




[(p) + (e)(m)]^u(t 1 z 1 



lim 



MNP l 



VA 



du 9va 



dz 



u - 



dz 



dt + uv Ao dz 



t dt J 



Si M et P coïncident on a (m) = (p), donc 
(m)[l + (e)]V7ru(ti*i) 



lim 



MNP 



va t, u ^r~ dt + uv a dz 

oz oz ' 



\L VAj ^n d<T 



Dans le second cas où (e) = —1 le premier membre est toujours nul. 
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Envisageons le second cas où (e) = 1. Alors 

du dvA 



2(m)-\/n u(ti z\) = lim 



MNP 



va ^ li- 
ez dz 



dt + uva dz 



-iL vAj ^ da > 



Supposons maintenant que la courbe MNP appartienne à un plan parallèle 
à xy, alors dz = sur la courbe, et l'on a 



2(m)\fTt u{tiZi) =\iml \ 

U mnp 



du dv A 
VA °-d- Z - U -dz- 



dt 



1 



f dt J 



Si / = la formule se réduit encore et devient 



2(m)y/7ru(ti z\) 



lim 



MNP 



VA 



du 
dz 



u - 



dv Ao 
dz 



dt 



Enfin si nous prenons pour MNP un cercle ayant le centre sur l'axe x la 
formule précédente devient la formule (3') de Poisson. 

De même les formules de Schlaefli peuvent se déduire comme cas particuliers 
de celles que nous avons données. 

Il faut remarquer que toutes ces formules se rattachent au premier cas, c'est 
à dire à celui où la projection de la surface a entoure le point xo,0. Or cela ne 
peut arriver que si l'on suppose que la courbe MNM' n'appartient pas au plan 
xz. Il fallait donc sortir des valeurs réelles de t pour trouver ces formules, et c'est 
justement ce que nous avons fait. Ce n'est qu'en sortant du domaine réel et en 
pénétrant dans le domaine complexe que la polydromie de la fonction (7) a pu 
jouer son rôle, car la polydromie ne ressort qu'en parcourant des cycles autour de 
la droite KK' . Nous avons mis en lumière qu'à cette polydromie se rattache tout 
le succès de la méthode nouvelle qu'on a employée. 

14. Nous finirons en donnant une application de la formule de Poisson, ce 
qui nous amènera à revenir sur un point auquel nous avons touché dans la 8 ieme 
leçon (§ 1). Nous avons dit qu'il y avait des problèmes de la théorie des ondes qui 
ne ressortent pas des équations du type hyperbolique. Ce sont les problèmes des 
ondes des liquides incompressibles. Si l'on envisage le cas le plus général des petites 
vibrations des liquides assujettis à des forces quelconques ayant un potentiel, on 
trouve que le problème se réduit à trouver une fonction V harmonique dans le 
domaine S occupé par le liquide qui le long de la surface libre u du liquide vérifie 
dV d 2 V 



l'équation 



-a 



dn 



vérifie l'équation 



dV_ 

dn 



dt 2 



et sur les parois rigides a qui renferment le liquide 



0, n étant la normale au contour dirigée vers l'intérieur 



de l'espace S, a une quantité positive, t le temps. 
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Le problème dépend de l'équation différentielle de Laplace 

A 2 V = 0, 

car V doit être harmonique. Il est donc rattaché à une équation du type elliptique. 
Ce sont les conditions le long de la partie u du contour qui déterminent le caractère 

ondulatoire du mouvement. 

fdV\ 
Or on démontre que V est déterminée si l'on connaît les valeurs Vn et — — = 

\dtj 

dV 
Vq de V et de — — sur u> pour t = 0. On peut aussi donner une méthode générale 

pour résoudre le problème. Soit U A (x,y,z,t) une fonction harmonique régulière 
de x,y,z dans le domaine S, excepté dans le point A de coordonnées xi,yi,zi où 
l'on a 

U A (x, y,z,t) = h W A (x, y, z, t), r A = y/(x 1 - x) 2 + (y x - y) 2 + (z x - z) 2 , 

r A 

W A (x, y, z, t) étant une fonction régulière. Si nous supposons que sur u soit vérifiée 
l'égalité 

dU A d 2 U A 



-a 



et sur a 



dn dt 2 

dU A 



dn 







dU A 



et pour t = t\ sur us U A et — - — soient nulles, alors on aura 

dt 



(12) Vixxyxzitx) 



1 d 
An dt 1 



(U. 



fdV 



A 



- 



V dt 



Va 



dt 



du, 



(\ 



(f)TT \ f)T T 

—^— I désignant les valeurs de U A et — — pour t = 0. Cette formule 

sert pour reconduire le problème de déterminer V à celui de trouver U A , c'est 
à dire à déterminer une fonction qui à certains points de vue est analogue à la 
fonction de Green. 

15. Mais laissons de côté la théorie générale qui nous conduirait trop loin, et 
envisageons le cas où le liquide occupe un domaine sphérique et où u> est toute la 
surface de la sphère. Supposons aussi que a soit constante de sorte que l'on puisse 

dV d 2 V 
prendre a = 1. Alors sur la surface de la sphère on aura 



dr 



dt 



2 ' 



r étant le 



76 



rayon vecteur conduit par le centre de la sphère. Posons logr = p, et supposons 
que le rayon de la sphère soit égal 1. Alors l'équation précédente pourra s'écrire 

%-%■ 

Cela posé, r, ■#, cp étant les coordonnées polaires soient V (r, fi, cp) et Vq(t, #, <p) 

les fonctions harmoniques à l'intérieur de la sphère qui sur la surface de la sphère 

dV 
deviennent égales aux valeurs V et Vq données de V et de — — pour t = sur 

la surface même. On pourra aussi les considérer comme des fonctions de p, $, <p, 
et l'on écrira Vo(p, #, <p), V^p,^, <p). Regardons ê et <p comme des quantités con- 
stantes, et déterminons par la formule de Poisson l'intégrale de l'équation (13) qui 
pour t = se réduit à Vo(p, $,cp) et dont la dérivée par rapport à t pour t = 
se réduit à V^p, i?, ip). Appelons cette intégrale V(p, #, <p, £). On voit facilement 
qu'elle est la fonction qu'on veut déterminer. 
En effet on aura 

(14) A 2 — = . 

1 ' dp dt 2 

Mais 

A 



2 dV 1 8(r 2 A 2 V) 



dp r 2 dp 
par suite l'équation (14) s'écrit 

8(r 2 A 2 V) r 2 8 2 (A 2 V) 



(15) 



dp dt 2 



d(r 2 A 2 V) 

Or r 2 A 2 V pour t = devient r 2 A 2 V = et — — pour t = devient 

dt 

t 2 A 2 Vq = : par suite r 2 A 2 V = 0. Donc V est harmonique, vérifie l'équation 

dV d 2 V dV 

— — — = -——- à la surface de la sphère et en même temps V et — — pour t = 

dn dt 2 dt 

prennent à la surface de la sphère les valeurs données, comme il fallait démontrer. 

Il est facile de vérifier que V est régulière au centre de la sphère. Le problème est 

donc résolu. 

On pourrait aussi se servir de cette méthode pour calculer la fonction U a 

analogue à la fonction de Green et calculer ensuite V par la formule (12). 
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du du 
Schlaefli. Uber die partielle Differentialgl. — ——, - = 0. Journ. de Crelle 

dt dx 2 
T. 72. 

Betti. Sopra la determinazione délie température nei corpi solidi omogenei. 

Memorie délie Soc. Italiana délia Scienza (detta dei XL). S. III, T. I. P° IL 

z dz 
Appell. Sur l'équation — - - — = 0. Journ. de Math. T. VIII, s. 4. 

dx 2 dy 

Volterra. Sur les équations différentielles du type parabolique. Comptes ren- 
dus des Séances de l'Ac. des Sciences. 5 Dec. 1904. 



Note additionnelle à la 5 ieme leçon. 

Communication faite le soir du 2 Mars. 
1. Rappelons la formule bien connue de trigonométrie 
sin(a + (3) = sin a cos (3 + sin (3 cos a, 
on peut en tirer que si 

(1) sin a cos/? + sin/? cos a = const. 

on doit avoir 

a + (3 = const. 

Posons maintenant 

f x dx [ y dy 



o Vl - x 2 Jo a/1 -y 2 



l'équation (1) s'écrira 

(2) x a/1 — y 2 + yVl — x 2 = const. 
Donc si cette égalité est vérifiée on doit avoir 

(3) P dX + r dy = const 

Jo \/l - x 2 Jo a/1 - y 2 



(4) 
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Rapportons-nous à deux axes coordonnés x,y, et soit s 
un arc de la courbe ayant pour équation (2), M un point de la 
courbe, M x et M 2 ses projections sur les axes ; alors OM x et 
OM 2 seront les amplitudes des deux intégrales qui paraissent 
dans la formule (3). 

Le théorème de trigonométrie peut s'énoncer en disant : 
si on déplace le point M sur la courbe s la somme des intégrales 




dx 



Vï 



x< 



dy 



y/ï^ 



étendues respectivement à OM 1 et OM 2 est constante. 

2. Passons à une extension de ce théorème au cas de trois variables. 
Envisageons les trois intégrales 



dydz 



A 



A 



3 2 

-y 2 



Pi 



dzdx 



-x" 



/./, 



6i 



z 2 + (3 2 



dx dy 



6i 



-îr 



-x z 



Ps 



où &i,&2,&3; fiiifhiflz; A, /i, // sont des quantités constantes et A + \x + v = 0. 
Prenons la surface algébrique ayant pour équation par rapport aux axes coordonnés 
x,y,z 



(5) \x 




6i 



z 2 + l3 2 + vz 



6i 



-x z 



fo 



const. 




et menons une ligne quelconque s fermée 
sur cette surface. Soient si, s 2 , s 3 les pro- 
jections de cette ligne sur les plans coor- 
donnés, et supposons que nous étendions 
les trois intégrales (4) aux aires a±, a 2 , a 3 
renfermées dans ces projections. Alors la 
somme des trois intégrales (4) ne chang- 
era pas en déplaçant et en déformant la 
ligne s sur la surface (5). 

La démonstration de ce théorème 
ne présente pas de difficultés. Il donne 
évidemment une extension du théorème 
élémentaire de trigonométrie par l'em- 
ploi des fonctions de lignes. La remarque 
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la plus intéressante à faire est que la surface algébrique constitue une liaison al- 
gébrique entre les trois lignes dont dépendent les trois intégrales. Ce théorème 
peut bien se généraliser. 



Voir Volterra : Un teorema sugli integrali multipli. Att. délia R. Accad. di 
Torino, 1897. 
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